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Then we prove that under natural assumptions, the solutions of the free boundary motion converge 
to solutions of the Euler equations in a domain with fixed boundary when the coefficient of surface 
tension tends to infinity. 

Contents 

1. Introduction 2 

1.1. Organization of the paper 6 

1.2. Scaling by length 7 

1.3. A geometric interpretation of theorem 1.2 7 

2. Auxiliary results 8 

3. The space S^{Ul) 11 

4. A new system of equations 12 

5. Geometry of the boundary and analysis of V/ 14 

5.1. A rewritten equation for V/ 15 

5.2. Local coordinates 16 

6. Analysis of f\ dn 19 

6.1. Analysis of / in the interior 37 

7. Proof of theorem 1.2: existence 38 

7.1. Overview of the argument 38 

7.2. Successive approximations 39 

7.3. Convergence of the approximating sequences 45 

7.4. Solution 47 

7.5. Proof of theorem 1.2: existence 48 

8. Proof of theorem 1.2: convergence 49 

8.1. Proof of corollary 1.8 53 

References 54 


E-mail addresses: marcelo.disconzi@vauderbilt.edu, ebin@math.sunysb.edu. 

Key words and phrases, large surface tension, free boundary, Euler equations, incompressible fluid. 
Marcelo M. Disconzi is partially supported by NSF grant 1305705. 

1 



2 


DISCONZI AND EBIN 


1. Introduction 

Consider the initial value problem for the motion of an incompressible inviscid fluid with free 
boundary whose equations of motion are given in Lagrangian coordinates by (see below for the 
equations in Eulerian coordinates) 


II 

1 

O 

in 17, 

(1.1a) 

div(it) = 0 

in 77(17), 

(1.1b) 

p\dr,m = kA 

on c?T7(17), 

(1.1c) 

. h( 0 ) = id, 77 ( 0 ) = u 0 , 


(lid) 


where 17 is a domain in R n ; r){t, •) is, for each t, a volume preserving embedding 77 (f) : 17 —>• R n 
representing the fluid motion, with f thought of as the time variable ( 77 ( 7 , x) is the position at 
time f of the fluid particle that at time zero was at x); denotes derivative with respect to f; 
17(f) = 77 (f) (17); u : 17(f) —> M 71 is a vector held on 17(f) defined by u = 77 o r 7 -1 (it represents the fluid 
velocity); A is the mean curvature of the boundary of the domain 17(f); p is a real valued function 
on 17(f) called the pressure; finally, k is a non-negative constant known as the coefficient of surface 
tension, id denotes the identity map, uq is a given divergence free vector held on 17, and div means 
divergence. The unknowns are the fluid motion 77 and the pressure p, but notice that the system 
( 1 . 1 ) is coupled in a non-trivial fashion in the sense that the other quantities appearing in ( 1 . 1 ), 
namely u, A, and 17(f), depend explicitly or implicitly on 77 and p. 

With suitable assumptions, we shall prove the following result, concerning the existence of solutions 
to ( 1 . 1 ) and the behavior of solutions when the coefficient of surface tension is large, i.e., in the limit 
k —>• 00 . A precise statement is given in theorem 1.2 below. 

Theorem (Main Result, see theorem 1.2 for precise statements). Under appropriate 
conditions on the initial condition uq and on <917, we have: 

1) If k> 0, then (1.1) is well posed. 

2) Consider a family of initial conditions uq k parametrized by the coefficient of surface tension 
that converges, when k —>• 00 , to a divergence free and tangent to the boundary vector field Dq. Then, 
the corresponding solutions to (1.1) converge to the solution of the incompressible Euler equations 
on the fixed domain 17, given by 

( = —Vtt o (, (1.2a) 

< div(CoC _1 ) = 0 , ( 1 . 2 b) 

. C(0) = id, C( 0 ) = T? 0 - ( 1 . 2 c) 

Here, Q(t, ■) is, for each t, a volume preserving diffeomorphism ((t) : 17 —>■ 17. 

Remark 1.1. It is well known that the pressure 7 r in the incompressible Euler equations is not an 
independent quantity, since it is completely determined by the velocity vector field 7? = ( o £ _1 (see, 
e.g., [26]). 

We remind the reader that in Eulerian coordinates, equations (1.1) and (1.2) take, respectively, 
the following forms: 

, ^+V n ti = -Vp in (J {f} x 17(f), 

0 <t<T 

< div(tt) = 0 in 17(f), 

p = kA on <917(f), 

M 0 ) = u 0 , 


(1.3) 
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and 


ff+ V*0 = -Vvr 
div(5) = 0 
(5, v) = 0 
^(0) = i?o, 


in [0, T] x 
in fl, 
on 30, 


(1.4) 


where in the free boundary case, p. k, A and vr, are as before, u is the velocity field with uq as its 
initial value. In the fixed boundary case u is replaced by 5 which has an initial value "do- The other 
symbols have the same meaning, except v which is the unit outer normal to 50. In this case, our 
theorem states that u converges to i? as k goes to infinity. 

In Eulerian coordinates, the free boundary Euler equations also carry a boundary condition stating 
that the normal speed of the moving boundary equals to ( u,N ), where N is the unit outer normal 
to 50(f). This is the same as saying that the vector field dt + V M is tangent to (J {f} x 0(f). 

o <t<T 

In order to state the main result, we need to introduce some definitions. Given manifolds M and 
N, denote by H S (M,N ) the space of maps of Sobolev class s between M and IV; that is, maps with 
derivatives up to order s in L 2 . For s > ^ + 2 define 

£®(0) = £* = 1 77 €E H s (yi,W l ) J(rj) = 1, 77 —1 exists and belongs to H s (r](Q), 0) j, 

where J is the Jacobian. £®(0) is therefore the space of IL s -volume-preserving embeddings of O into 
R n . Define also 

XT)(0) = D* = e H s { 0,R") J{rf) = 1,7/ : O —y O is bijective and j? - 1 belongs to H s |, 
so that D^(D) is the space of i5 s -volume-preserving diffeomorphisms of D. Notice that C 

W- 

Let B s 5 + 2 (d£i) be the open ball about zero of radius Jo inside H s+2 (d£l). We shall prove that if Jo 
is sufficiently small, then the map 


where / satisfies 


v ■■ #! 0 +2 (^) 
<p{h) = f, 

J(id+V/) = 1 
f = h 


H s+ I(n), 


in D, (1.5a) 

on 5D, (1.5b) 

is a well defined C 1 map, and is a smooth submanifold of H s+ 2 (D). 

We note that the map (p solves a non-linear analog of the Dirichlet problem, that of extending h 
from 5D to a function on D. In fact if D C R 3 with standard coordinates then (1.5a) can be written 

A/ + fxxfyy + fxxfzz + fyyfzz ~ fly ~ fl z ~ fy Z + det (D 2 f) = 0, 

so the difference between (1.5) and the Dirichlet problem are the non-linear terms, which will be 
shown to be small. The purpose of (1.5a) is to ensure that id +V/ is volume preserving. 

Using tp we then construct another map 

$:P^(D)x^|+ 2 (5D))^^(D), 
defined by <h(/3, /) = (id +V/) o (5. 

Thus 4>(/3, /) is the composition of two volume preserving maps. 


( 1 . 6 ) 
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We define C £^(£i) by 

= <${v;(n)xv(B s + 2 (dn))). 

Notice that since j3 € D^(Q), we have j3(d£l) = dQ. Therefore, solutions rj to (1.1) that belong to 
decompose to a part fixing the boundary and a boundary displacement, i.e., 

r/ = (id +V/) o f3. (1.7) 

The decomposition (1.7) is one of the main ingredients of our proof, and establishing that rj —>• £ as 
k —} oo will be done by showing that the boundary displacement V/ goes to zero as k —>• oo. We 
shall also show that, under our hypotheses, V/ is in fact | degree smoother than rj (though /3 is as 
regular as 77 ). In this sense, we work with embeddings which have smoother boundary values. A 
more detailed discussion for the motivation for introducing S'MM) is given in [19]. 

We are now ready to state our main result. The usual decomposition of a vector held X into its 
gradient part, QX, and divergence free and tangent to the boundary part, PX, which appears in 
the hypotheses of the theorem, is reviewed in section 2. We note that if X is not tangent to the 
boundary, we will have QX 0 even if div(X) = 0. We let || ■ || s denote the Sobolev norm. We 
shall state and prove the theorem only in three-spatial dimensions, since this is the case of primary 
interest. We point out, however, that basically the same proof works in higher dimensions, except 
that the calculations increase significantly in complexity. 

Theorem 1.2. Let s > 4 + 2, Ll be a bounded domain in R 3 with a smooth boundary and uq € 
H S (D,M. 3 ) be a divergence free vector field. Denote by Quq the gradient part ofuo. 

1) If k > 0, then there exist a T K > 0 and a unique solution (jg K ,p K ) to (1-1), with initial condition 
uq. The solution satisfies: 

Vk € C°([0,T k ),^(O)), r, K € L°°([0,T K ),H s (M)), t ) k € L°°([0, T K ), H s ~l (H)), 
p K € L°°([ 0 ,T K ),i 7 s_ 2 (n K (t))), where VL K {t) = r] K (t)(Q). 

Moreover, for each t € [0, T K ), 

g K {t) is in <^(fi). 

2) Let {uok} C H s (n,R 3 ) be a family of divergence free vector fields parametrized by the coefficient 
of surface tension k, satisfying || Quo K || s < for some constant C, and such that uo K converges 
in H S (Q, R 3 ), as k —>• 00 , to a divergence free vector field i9q which is tangent to the boundary. Let 
C € C' 1 ([0,T],'D*(H)) be the solution to (1.2) with initial condition Dq, defined on some 1 time interval 
[0,T]. Assume that the mean curvature of dQ is constant, and let (■ r] K ,p K ) be the solution to (1.1) 
with initial condition uq k and defined on a time interval [0, T K ), as stated in part (1) above. Finally, 
assume that [0, T K ) is taken as the maximal interval of existence for the solution ( t!k,Pk )• Then, ifT 
is sufficiently small, we find that T K >T for all k sufficiently large, and rj K (t) —>• C(f) as a continuous 
curve in £*(fl) as n —>• 00 . Also, rj K (t ) —> C(f) in H S (D) as k —» 00 . 

Remark 1.3. We stress that in the theorem exists and is unique by [26]. 

Remark 1.4. Note that since rj K (t) € Soffit), in particular g K satisfies decomposition (1.7), and £l(t) 
has a i7 s+1 -regular boundary. 

^Notice that [0,T) is not a maximal interval of existence since the solution ( exists on the closed interval [0,T]. It 
can however be arbitrarily close to the maximum. 
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Since existence of solutions to (1.1) has already been established in the literature (see [10, 40] 
and comments below), the main result of this paper is the study of the singular limit k —>• oo and 
the corresponding convergence of solutions, i.e., part (2) of theorem 1.2. It should be stressed that, 
in this regard, the hypothesis of constant mean curvature at time zero cannot he removed. Indeed, 
on physical grounds, if the mean curvature Aqq of 3D is not constant, one expects that f l(t) will 
develop high-frequency oscillation for large k, but solutions will not converge in the limit k —>• oo. 
This is because the dynamics can be thought of as mimicking the behavior of motions with a 
strong constraining force, as we explain in section 1.3 (see also remark 6.1). We also notice that, 
from the point of view of k —>• oo, the assumption that Quo be small in part (2) of theorem 1.2 is 
natural, since Quo has to be small if uq is near ido- 

Despite the fact that, as mentioned above, well-posedness of (1.1) is known in the literature, our 
methods are entirely different than previous works, thus of independent interest. 

Remark 1.5. In this paper we treat exclusively the three-dimensional case, but the same proof 
works in n = 2. In fact, the calculations of section 5 simplify in two dimensions. 

Remark 1.6. In terms of the more familiar Eulerian coordinates, theorem 1.2 asserts, in particular, 
the convergence of u K (t) o r) K (t) to id(t) o £(t) in H S (Q). We believe, however, that in this problem 
the statement in Lagrangian coordinates looks more natural. This is because u K and id are defined in 
different domains, and therefore only the convergence of u K o ij K to id o £, and not of u K to id, makes 
sense. In other words, to state the convergence in Eulerian coordinates, we still need to invoke the 
flows and (. 

Remark 1.7. We point out that, due to a classical result of Aleksandrov [2], if Aqq is constant 
then 3D is a sphere. This fact, however, is not used directly in our proof. Had we not known of 
Aleksandrov’s result, our proof would still follow solely from the constancy of Adu- In particular, we 
expect theorem 1.2 to hold in the more general situation where 3D is the boundary of a region inside 
a Riemannian manifold with constant mean curvature with respect to the corresponding Riemannian 
metric. 

Theorem 1.2 gives some interesting insight into the structure of solutions to the free boundary Euler 
equations: by uniqueness, any solution to (1.1) satisfying our hypotheses will obey decomposition 
(1.7), regardless of the method employed to construct such solutions. 

As already pointed out, one wishes to show that V/ K in decomposition (1.7) goes to zero when 
k —^ oo in order to establish theorem 1.2. Since in (1.7) j3 K € D^(D), one expects that not only does 
—>• (, but fi K —>• C as well. This is in fact the case: 

Corollary 1.8. With the same assumptions and notation of theorem 1.2, consider the convergence 
r] K C- Then we also get /3 K (t) f(t) and /3 K (t ) —> ((t) in H s ( D), V/ K (t) —>• 0 in H s+ 2 (D), and 
Vf K {t) 0 in H s ( D). 

It is interesting to note that while r] K —>• ( and rj K —>• f, in general the corresponding pressures do 
not converge, even if the initial data are C°°. To see this, we first point out that since n is defined 
up to an additive constant, and thus only V7r is well-defined, one can only speak of convergence of 
Vp K to V7r, and not of p K to 7 r. Consider the two-dimensional case for simplicity, pick any function 
/ which is constant on 3D and let uq = (f y ,—f x )- Then uo will be divergence free and tangent to 
the boundary. The pressure for (1.2) at time zero will then satisfy: 

— An = 2(f xy — f X xfyy)i 

and V i,n will equal zero on 3D. Thus 7r in general will not be constant on 3D, so one cannot expect 
that Vp K , the solution of (1.1) will converge to V7r, as k —>• 00 , even at time zero. As a consequence, 
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convergence of the second time derivatives, i.e., rj K —>• C, generally fails (see the analogous results 
in [23, 25]). However, due to the convergence of the first time derivatives, the pressures over any 
positive time interval converge: 

Corollary 1.9. Under the same assumptions and notation of theorem 1.2, one has ff b Vp K o rj K —>• 
f* b V7T o ( in H S (Q), for any 0 < t a < 4 < T. 

The convergence part of theorem 1.2, namely, part (2), was proven by the authors in two spatial 
dimensions in [19]. However, compared to [19], theorem 1.2 is self-contained, in the sense that the 
existence of rj K is established before proving the convergence rj K —> (, whereas in [19] we relied on 
the existence results in Coutand and Shkoller [10] in order to obtain the convergence. 

The mathematical study of equations (1.1) has a long history, although for a long time results 
only under restrictive conditions had been achieved. In particular, a great deal of work has been 
devoted to irrotational flows, in which case the free boundary Euler equations reduce to the well- 
known water-wave equations. See [3, 4, 8, 39, 44, 66, 68]. More recent results addressing the question 
of global existence can be found in [30, 32, 67] and references therein. 

Not surprisingly, when equations (1.1) are considered in full generality, well-posedness becomes a 
yet more delicate issue, and most of the results are quite recent. In this regard, Ebin has showed that 
the problem is ill-posed if k = 0 [21], although Lindblad proved well-posedness for n = 0 when the so- 
called “Taylor sign condition” holds [40, 42]; see also [9] (the linearized problem was also investigated 
by Lindblad in [41]). When n > 0, a priori estimates have been obtained by Shatah and Zeng [53], 
with well-posedness being finally established by Coutand and Shkoller [10, 11] (see also [49]). See 
also [54]. Coutand and Shkoller also established the convergence of solutions in the limit k —>• 0 + 
in the case that the Taylor sign condition holds. Other recent results, including the study of the 
compressible free boundary Euler equations and singularity formation, are [6, 7, 12, 13, 14, 15, 16, 31]. 
We point out that the analogous free boundary problem for viscous fluids was first and extensively 
studied by Solonnikov [43, 55, 56, 57, 58, 59, 60, 61], with some more recent advances found in 
[37, 47, 50, 51, 52] and references therein. 

Lindblad’s result [40] is based on a Nash-Moser iteration, while Coutand and Shkoller [10] obtained 
existence and convergence when k —>• 0 + by developing a technique they call convolution by layers. 
The results here presented, based on the decomposition defined by (1.6), provide yet a third, different 
method of proof (valid for k > 0). We point out that, with exception of the authors’ work [19] in 
two-dimensions, the limit k —>• oo does not seem to have been investigated in the literature before. 

Notation 1.10. We reserve 12 for the fixed domain, with 12(f) being always the domain at time t, 
i.e., 12(f) = r](t)(Q). Of course, 12(0) = 12. In several parts of the paper the subscript k will be 
dropped for the sake of notational simplicity. 

Notation 1.11. H s {i 2) and iL s (<9I2) denote, respectively, the Sobolev spaces of functions on 12 and 
<912, with norms || ■ || s and || • || s g. LP(I2,M n ) etc are similarly understood, although when the 
manifolds are clear from the context, we simply write H s or H S (M ) for H S (M, N). Notice that 
H° denotes the L 2 space, with norm || • ||o- H q denotes the Sobolev space modulo constants. We 
use both V and D to denote the derivative. D w is the directional derivative in the direction of 
w, w a vector. The letter C will be used to denote several different constants that appear in the 
estimates. Sometimes we write C = C (a, b,...) to indicate the dependence of C on a, b ,.... We use 
the following abridged notation for partial derivatives: -£j = <9,;, & 'f (1 -, = <9 = dijk, etc. 

1.1. Organization of the paper. This paper is organized as follows. In section 1.2 we make some 
remarks about the role of surface tension. In section 1.3 we give a geometric interpretation of our 
theorem in terms of curves in the group of volume preserving diffeomorphisms and volume preserving 
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embeddings. In section 2, we state several known results that will be used and fix some notation. 
In section 3, we carry out the construction of the space <S)f (12) as outlined in the introduction. In 
section 4, we derive a new set of equations that splits the dynamics into an equation for a function / 
that controls the boundary motion, and an equation for a diffeomorphism f3 that fixes the boundary 
setwise. / is determined by its boundary values. Thus in section 5, we derive a further equation for 
f\on- Section 5 consists solely of a series of calculations necessary to analyze f\ a ^ and some readers 
may want to skip it. Section 6 is the core of the paper, where the existence of f\ afl is established 
and estimates for f in terms of - are obtained. Section 7 establishes the existence of solutions to 
(1.1) via an iteration scheme. Section 8 shows the convergence t/ k —>• £ in the limit k —>• oo. 

Remark 1.12. Throughout sections 3 to 7.4 we work under the hypotheses of part (2) of theorem 

1.2, i.e., we assume that k is large, <912 has constant mean curvature, and || Quo || s < In section 
7.5 we show how to prove the existence result under the general assumptions of part (1) of theorem 

1 . 2 . 

1.2. Scaling by length. When we speak of large surface tension or large k, we should take into 
account the size of the domain 12. It seems clear that surface tension should have more of an effect 
in a small domain than in a large one. To clarify this we examine the effect of scaling the length of 
the domain. 

Let A be a positive scale factor and assume 7/(2) is some motion satisfying (1.1). Then on the 
scaled domain A12 define C(2) by £(2)(Ax) = \rj(t)(x). Then letting y = Xx we find that 

C(2)(y) = A 7/(2) (x). 

A routine computation shows that ( satisfies (1.1) on A12 with p replaced by q where q is defined 
by q{y) = A 2 p(x). However the mean curvature of <9C(A12) = <9A?/(12) is (1/X)A 1 where A is the 
mean curvature of <97/(12). Thus q = X 2 p = X 2 uA = A 3 k( 1/A)AI so the scaled motion has an effective 
coefficient of surface tension of A 3 k. Hence, when we study the effect of surface tension we really 
should consider u divided by a typical length cubed or k divided by the volume of the domain. For 
a given k, the surface tension will have a much greater effect on a drop of liquid than it will on a 
large body. 

1.3. A geometric interpretation of theorem 1.2. Theorem 1.2 not only gives a satisfactory 
answer to the natural question of the dependence of solutions on the parameter /c; it also addresses a 
well motivated problem in Applied Science, namely, when one can, by considering a sufficiently high 
surface tension, neglect the motion of the boundary in favor of the simpler description in terms of 
the equations within a fixed domain. 

The physical intuition behind theorem 1.2 is very simple, as we now explain. The system (1.1) 
can be derived from an action principle with Lagrangian 

£( 7 /) = K{rj) - H(r/), (1.8) 

where 

K (V,V) = \ \h? (1-9) 

is the kinetic energy and 

V(rj) = tc|<912(2)| — ac|cA 2(0)| = Area(<912(2)) — Area(<912(0))^ 


( 1 . 10 ) 
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is the potential energy 2 , and with rj : [0, T) —> £*(17). The energy for the fluid motion (1.1) is given 
by the sum of the kinetic and potential energies (1.9) and (1.10), respectively, 

E(t) = K(i 1 ,r l ) + V( V ) 

= f a \v\ 2 + k (\dv(t)\-\dn i), ( ' 1 ' 11 ) 

This energy is conserved, and therefore 

m = \ ii «o iio (i-i2) 

where we have used r) = u o g and rj(0) = id. 

Our theorem 1.2 is almost an example of a general theorem on motion with a strong constraining 
force [23]. For the general theorem we are given a Riemannian manifold M and a submanifold N. 
Also given is a function V : M —> M. which has N as a strict local minimum in the sense that W = 0 
on N and D 2 V is a positive definite bilinear form on the normal bundle of N in M. Then if g K (t) is a 
motion given by the Lagrangian C(g, rj) = \ ( 77 , 77 ) — KV(g) where (, ) is the Riemannian metric, and 
if £(f) is a Lagrangian motion in N of ^(C,C) with the same initial conditions as g K {t ), the theorem 
says that g K {t) converges to ((t) as k —>• 00 . Also rjk ~> (, but the second derivative in general does 
not converge. For our theorem, M = £*(!7), N = P*(17), ( , ) is the L 2 inner product on tangent 
vectors and V(rj) is given by (1.10). 

Our theorem 1.2 is not actually an example of the general theorem for two reasons: 

a) The L 2 inner product on tangent spaces is only a weak Riemannian metric 3 . The topology that 
it induces is weaker than the H s topology of £*(17). 

b) The bi-linear form D 2 V is only weakly positive definite on each normal space; it gives a topology 
weaker than the H s topology. 

Thus, while theorem 1.2 is not a particular case of established results about the behavior of the 
Euler-Lagrange equations near a submanifold which minimizes the potential energy [23], it can be 
viewed to be in the spirit of those results. Here, as in [23, 24, 25], the manifold minimizing the 
potential energy is P*(17); see also [18]. 

2. Auxiliary results 

Here we recall some well known facts which will be used throughout the paper. For their proof, 
see e.g. [1, 5, 17, 22, 46]. 

Proposition 2.1. Let s > + 2, g E f € H s (g(Ll)). Then f o g £ H s (Ll) and 

II /09 !U< c II / II, ( 1 + II g II,), (2.1) 

where C = C(n, s, £2). 

We shall make use of the following well-known bilinear inequality 

|| uv || r < C\\u || r || v || s , ( 2 . 2 ) 

for s > s > r > 0, where C = C(n,s,r,fl). 

2 Many authors consider instead V(rj) = «;|9fi(£)|. As the equations of motion remain unchanged by adding a 
constant, we choose to normalize the potential energy to make V = 0 at time zero. Such a normalization is convenient 
for our purposes as we are interested in taking k —¥ 00 , in which case, if we did not subtract the contribution at time 
zero, V ( 9 ) would diverge to infinity. 

3 We recall that a weak Riemannian metric is one which induces, on each tangent space, a weaker topology than the 
one given by the local charts. This is a feature exclusive to infinite dimensional manifolds; see [20, 22, 38] for details. 
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For future reference, we remark that (2.2) still holds true in negative norm Sobolev spaces (in a 
compact domain without boundary). Indeed, if a € H~ ri , n > 0, b e FT' 2 , r 2 > §, then 

|(ab, oj)q| _ |(a,Mo| 


ab |_ ri = sup 
< sup 

u&H r l 


OJ 


Ul 


a 


r l 


= sup 
0j£H r l 

boj 


OJ 




|ri 


< c 


a _ 


OJ 


r i I 


in > 


in 


after using (2.2) in the last step to estimate || boj || ri < C || oj 

Notation 2.2. Although the dimension n = 3 is fixed throughout, we sometimes write n instead of 
3 in order to make it easier to read off conditions such as s > | that are necessary for the application 
of (2.2) and other dimension dependent results. 


Recall also that restriction to the boundary gives rise to a bounded linear map, 


TI ||s,<9— C 


ls+ 


i, s > 0, 


(2.3) 


with C = C(n,s,Q). Estimates (2.2) and (2.3) will be used throughout the paper, so we shall not 
explicitly refer to them in every instance. 

The following div — curl estimate is well-known (see, e.g., [63]): let II be a domain with an H r 
boundary, r > 3, v be a vector field on 12 such that v € H°( 12), curlu € FF S_1 (I2), divu € FP _1 (I2), 
and (v, v) G H s ~ 3 (<9f2), where v is the unit vector normal to <912. Then, v € H S (Q), and we have the 
following estimate 

II v || s < (7(11 v ||o + || curl v |U_i + || div v || s _i + || (v,v) || s _i s ). (2.4) 

Next we recall the decomposition of a vector field into its gradient and divergence free part. Given 
an H s vector held oj on 12, define the operator Q : H s (i 2,M n ) —> ViF s+1 (f2,M n ) by Q{oj) = \7g, 
where g solves 

/A g = div(w) in f2, 

IS = on dn - 

Since solutions to the Neumann problem are unique up to additive constants, X7g is uniquely de¬ 
termined by oj, so Q is well defined. Define P : H s (i 2,M n ) —>• div _1 (0) I/ , where div _1 (0) y denotes 
divergence free vector fields tangent to <9f2, by P = I — Q, where I is the identity map. Then Q and 
P are orthogonal projections in L 2 . 

We shall make use of the following: 

Notation 2.3. If ?y : 12 —>• M n is a sufficiently regular embedding and O is a pseudo-differential 
operator defined on functions on r]{i 2), we let O v , which acts on functions defined on 12, be given by 

O v {h) = (0(h o r/ -1 )) or/. 


We remark that using notation 2.3, equations (1.1) can be written as 


^3: 

II 

1 

< 

o 

-3 

in 12 , 

( 2 . 6 a) 

div r? (? 7 ) = 0 

in 12 , 

( 2 . 6 b) 

Q <90 = 

on dQ, 

( 2 . 6 c) 

. 7 ( 0 ) = id, 77 ( 0 ) = u 0 , 


( 2 . 6 d) 


where q = p og and B = A op. Equations (2.6) reveal yet another advantage of Lagrangian 
coordinates, as all equations are now written in terms of the fixed domain 12. 
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Finally, theorems related to one-parameter groups of operators and abstract differential equations 
will be needed. We also state them here for the reader’s convenience. 

Theorem 2.4. Let X be a Banach space and let Z be a densely defined closed operator on X. 
Assume that every real A is in the resolvent of Z and that 

Hz + xr\\<fi y 

for some constant c\ > 0. Then, Z generates a C° semi-group of transformations e tz : X —> X, 
such that || e iz ||< ci. If Z' is a bounded operator with norm || Z' ||< C 2 , then Z + Z' also generates 
a C° semi-group and 

II e tZ ||< cie C2 ^, 


Theorem 2.5. Let X and Y be Hilbert spaces such that Y is densely and continuously embedded in 
X. Let {Z(t) | 0 < t < T} be a family of operators on X, each of which generates a C°-semi-group 
and assume that: 

(i) There exist constants a and f3 such that, for each t and for all positive t, 

|| e rZ(t) ||< ae 0T , 

(ii) For each t, e rZ restricted toY is a C° semi-group on Y. There exist constants 7 and 5 such 
that, for each t, there exists an inner product on Y, whose norm t || • || gives the topology ofY, and 
such that 

II e rZ(t) ||op( t) < Se 1T , 

where Op(f) is the operator norm on B(Y ) induced by t \\ ■ ||. Furthermore, there exist constants p 
and v such that, for any t\,t 2 € [0, T] and any y € Y, 

tall y II < ve ult2 ~ hl tl \\ v II- 

(in) Y is included in the domains of each Z(t), and Z{t) is continuous as a map from [0, T] to 
B(Y, X). 

(iv) Z(t) is reversible in the sense that Z(t) = —Z{T — t) also satisfies (i), (ii), and (Hi). 

Then, there exists a unique family of operators U(t,r ) € B(X), defined for t,r € [0, T], such that 

(a) U(t,r) is strongly continuous as a function of r and t, U(t,t ) = / (the identity operator), and 

II U(t,r) ||< cie C2|f - r| , 

for some constants ci,C 2 depending only on a, (3, 7 , 5, p, and u. 

(b) U{t,r) = U{t,T')U{r', t). 

(c) For all y €Y, 

( U(t,r)y ) = Z(t)U(t,r)y, 

where means right derivative at t = 0 and left derivative at t = T. 

(d) For all y 6 f, 

( U(t,r)y) = —U (t, r )Z (r )y, 

where means right derivative at t = 0 and left derivative at t = T. 
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(e) U(t,r)Y C.Y and for any t, rG [0,T], 

II U(t,T) ||o p (t)< c 3 e C4T e C5|t ~ r| , 

for some constants 03 , 04,05 depending only on a, f3, 7 , <5, fi, and v. 

(f) U(t,r) is strongly continuous into Y, as a function oft and t, and therefore 

^( u (t, T )y) = z (t)U{t,r)y 
is continuous in X as a function of t and r. 

Theorem 2.6. Let X, Y and Z(t) be as in theorem 2.5. Let eft) be a continuous curve in Y and 
define yft) to be 

y(t) = U(t,0)y 0 + [ U{t,T)e{r)dT, 

Jo 

for 0<t<T,yo£Y. U is the evolution operator given by theorem 2.5. Then yft ) is a continuous 
curve in Y which is C l in X, and it is the unique solution to the equation 

y(t ) = Z(t)y(t ) + eft) 

such that y( 0) = 7/0 - 

The proof of theorem 2.4 can be found, for instance, in chapter 12 of [29]. Theorems 2.5 and 2.6 
are proven in [33] (see also [34]). We remark that the results of [33] are much more general than the 
above. Here, we stated them in a form suitable for our purposes. See also [23]. 

3. The space S’f(Q) 

Here we construct the space S‘f 1 {iT), as outlined in the introduction. We assume that s > % + 2. 
To start we note that the equation 

J(id+V/) = 1 

can be written as 

Af+N(f) = 0, (3.1) 

where 

A/(/) = fxxfyy + fxxfzz + fyyfzz ~ fly ~ fl z ~ fy Z + det(T> 2 /). (3.2) 

Equation (3.1) can be considered as a non-linear Dirichlet problem for /, and so for / small, / should 
be determined by its boundary values. We shall present the argument for three dimensions, which 
is the main case of interest in this work. The interested reader can generalize the construction of 
SffLl) to higher dimensions. 

Given h € H s+2 (dLt), we are interested in solving 

fA/+AA(/) = 0 in n, 

\ f = h on dQ. 

F : H s+2 {dLl) x H s+ §(fi) -a H s+2 {dtl) x H s+ ^(Ll), 
by F(h,f) = (f\ dn -h,Af+M(f)). 


Define a map 


(3.3a) 

(3.3b) 
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Notice that F is C 1 in the neighborhood of the origin and F(0, 0) = 0, where we denote by 0 the 
origin in the product Hilbert space H s+2 (dLl) x H s+ 5 (17). Letting w € H S+2 (Q), we obtain 

D 2 F(0, 0)(w) = (w | 9n , Aw), (3.4) 

where D 2 is the partial derivative of F with respect to its second argument. From the uniqueness 
of solutions to the Dirichlet problem it follows that D 2 F(0, 0) is an isomorphism, and therefore by 
the implicit function theorem there exists a neighborhood of zero in H s+2 (dLl), (which we can take 
without loss of generality to be a ball B^ 2 (dfl)) and a C 1 map ip : B^ 2 (dfl) —> H s+ 5 (H) satisfying 
<£>(0) = 0, and such that F(h, pfih)) = 0 for all h € H^J" 2 (H). In other words, / = pfih) solves (3.3). 

Furthermore, Dtp = —(D 2 F)~ 1 DiF. Thus Dip is injective at the origin (in fact, it is not difficult 
to see that the derivative Dip( 0) is the harmonic extension map), and so ip is injective near zero. 
From this and the above it then follows that p(B s ^ 2 [dD )) is a submanifold of H s+ 2 (H). 

Recall now the definition (1.6). Notice that 4> is well defined (if ho is small) and its image belongs 
to £^(fl) since J(/3) = 1 and, by construction, J(id+V/) = 1. 

We have therefore proven: 

Proposition 3.1. Let s > + 2 and let B^ 2 (dQ) be the open ball of radius ho i n H s+2 (dLl). Then, 

if is sufficiently small, there exists an embedding ip : B'D 2 (dil) —> H s+ 2 (Ll), given explicitly by 
<p(h) = f, where f solves (3.3). Moreover, the map <J> given by (1.6) is well defined. 

Definition 3.2. Under the hypotheses of proposition 3.1, we define <^(H) C £^(fl) by 

^(n) = $(v( l (n)xip(B^ 2 (dn))). 


4. A NEW SYSTEM OF EQUATIONS 

In this section, we shall derive a different set of equations for the free boundary problem (1.1). In 
what follows, we shall make use of the well known decomposition of a vector field into its gradient and 
divergence free parts, as presented in section 2. Hence, recall that Q : H s (Q,R n ) —> V1L S+1 (H, R n ) 
and P : H s (Ll,R n ) —>• div^ 1 (0) y (where div _ 1 (0)^ denotes divergence free vector fields tangent to dLl) 
are the operators realizing this decomposition. They satisfy P + Q = I, where I is the identity map, 
and since VFP +1 (fl, M ra ) and div _ 1 (0) y are L 2 -orthogonal, it follows that Q and P are orthogonal 
projections in L 2 . 

To derive the new system, assume that solutions 77 to (1.1) can be written as 

r) = (id +V/) o (4.1) 

with (3 € V/ € H S (Q) and with / satisfying (3.3a). In this case we also observe that 

m = id, v/(0) = 0. 

It is customary to write the pressure as a sum of an interior and a boundary term, namely, 


P = Po + kA h 


( 4 . 2 ) 
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so that the system ( 1 . 1 ) takes the form 




' 77 = —Xp 0 p = — (Xpo + kXAh) 0 V 

in D, 

(4.3a) 


div(?) 0 r]~ l ) = 0 

in 77 (D), 

(4.3b) 


Ap 0 = -di v(X u u) 

in 77 (D) , 

(4.3c) 

< 

Po\d v (Q) = 0 

on 577 (D), 

(4.3d) 


A Ah = 0 

in 77 (D), 

(4.3e) 


Ah = A 

on 577 (D), 

(4.3f) 


„ 77(0) = id, j)(0) = u 0 . 


(4.3g) 

Differentiating (4.1) in time gives 




fj = (X f + v ■ DXf + v) 0 /3, 


(4.4) 

where v is defined by 

f3 = v 0 (3. 


(4.5) 

Using 77 ( 0 ) = id and 77 ( 0 ) = uq, from (4.4) we obtain 




uo = V/( 0) +n 0 . 



We therefore obtain 

O 

II 

h tl 

£ 

O 



and 

V/( 0) = Qu 0 . 



Differentiating (4.4) again and using (4.3a) gives the following equation: 


V/ + 2D v Xf + D 2 vv Xf + (v + v Xv)D Xf + v + v 

Xv = —Xpo (id +V/), 

(4.6) 

where the operator D' 2 v , 

acting on a vector w, is given in coordinates by 



(■ D 2 v w ) l = v J v l djdiw\ 


(4.7) 

or in invariant form by 

D 2 v w = D v X v w - D^ vV w. 



Define L on the space of maps from D to M n by 




L = id+D 2 /, 


(4.8) 

and let 

II 

fa 

Jrs 


(4.9) 

and 

fc-< 

to 

II 

0 

Jrs 


(4.10) 


where P and Q are as in section 2. Notice that L\ is invertible on the image of P if / is small, since 
in this case it will be close to the identity. Then we can write an arbitrary vector field X as 

X = LL^P(X) + {Q- L 2 L^P)(X), 

and thus effect the decomposition 

H s (Q,R n ) = LP{H s (Q,R n )) © Q(H s (n,R n )). 


(4.11) 
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Decomposing (4.6) in this fashion, and recalling that / also has to satisfy (3.1), we obtain 
' V/ + (Q- L 2 L^P)(2D v Vf + D 2 vv Vf + LQ(V v v )) 

= -(Q-L 2 L^ 1 P)(Vpo (id+V/)) in D, (4.12a) 

A/+AA(/) =0 in D, (4.12b) 

v + P{V v v) + L^P(2D v Vf + D 2 vv Vf) + L^P(LQ(V v v)) 

= -L^ 1 P(Vpo(id+V/)) in D, (4.12c) 

, V/(0) = 0, V/(0) = Qu 0 , v(0) = Puo, (4.12d) 

where M is given by (3.2). Equation (4.12c) implies that v also satisfies div(u) = 0 and ( v,u) = 0, 
where v is the unit outer normal to 3D, as v is in the image of P. The system (4.12) contains two 
equations for /, namely, (4.12a) and (4.12b), but no boundary condition. We shall transform these 
into a more standard evolution equation for / by restricting (4.12a) to 3D, obtaining an evolution 
equation for /| af2 , with an extension to D given via (4.12b). 

It is illustrative to point out that (4.12c) formally reduces to the Euler equations in the fixed 
domain when V/ = 0 and thus, formally, ( = /3. This is in agreement with the intuition discussed in 
the introduction that if the boundary displacements controlled by V/ approach zero when k —>• oo, 
the solution r] should approach £. Such an intuitive appeal notwithstanding, equation (4.12c) will not 
be used directly in our proof. There are two reasons for this. First, as explained in the introduction, 
there is no reason to suspect that the pressure will converge when we take the limit n —>• oo, hence 
no good control of the right-hand side of (4.12c) can be expected. Second, in view of the regularity 

3 

of p stated in theorem (1.2), the right hand side of (4.12c) will be in H s ~ 2 . Known techniques, 

3 

therefore, can only yield v, and hence (3, in H s ~ 2 . However, /3 € H s (£l) is required for 7 / € LP(D), 
see (4.1). Our proof overcomes these difficulties by making the most of the Lagrangian description of 
the fluid, which is consistent with the general idea that Lagrangian coordinates are “better behaved” 
than Eulerian ones (see similar discussion in [23] and [25]). 

5. Geometry of the boundary and analysis of V/ 

In light of proposition 3.1, f K is determined by its boundary values, provided it is small. In this 
section, we shall show that f K \gn obeys an equation of the form 

f fn — , p K , fQ T &k((3ki ^kiPk,i Ik) T (/^kj Vk • Pk ) On 3D 

\/*(0) = 0, f K (0) = h, 

where si K is a third order pseudo-differential operator on f K , is first order, ( € K is a zeroth order 
operator on v K , and fo and fi are known functions. The desired equation will be equation (6.1) 
below, and the present section is a derivation of (6.1) from (4.12a). This amounts essentially to a 
series of lengthy calculations, and some readers may want to skip them and move directly to section 
6 . 

From now on, the subscript k will be omitted in all quantities. Throughout this sections we 
assume we are given a sufficiently regular solution of (4.12). Moreover, as many of the derivations 
below are valid provided that / is sufficiently small, we shall assume so throughout this section. This 
smallness condition will be made precise in section 6. Several standard geometric constructions will 
be employed below. They can be found, for example, in [28, 62], 

Notation 5.1. Let 

rj = 77 o /fG 1 = id +V/. 
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Unless stated otherwise, from now on quantities with are defined on the domain r](£l). For example, 
if N denotes the normal to drj(£l), then N is the normal to drj(Q). 

5.1. A rewritten equation for V/. From (4.8), (4.9), and (4.10), we see that 

Q ~ L 2 L^P = Q (id - (id +D 2 f)L^P) 

= Q- QD 2 fL~ l P , 

where we used the fact that Q vanishes on the image of P. 

Let 

F = Ah ° V- (5-1) 

and 

qo=Po° V- ( 5 - 2 ) 

Notice that F depends on /. A calculation gives 

VAh °V = VF(Drf)~ 1 , (5.3) 

and 

Vpoorj = Vq 0 (Drf)~ 1 , (5.4) 

Then (5.1), (5.2), (5.3), (5.4), and (4.2) give 

Vp o rj = Vgo + V<7o ((Hi/) -1 — id) + kVF + kWF(— id +(Drf)~ 1 ). 

Using the above, (4.12a) can be written as 

V/ + kVF + kVA” 1 div (VF(- id +(Dr})~ 1 )^ - kVA" 1 div [d 2 fL^ l P{yF{Dy)~ l ) s j 
+ 2VA” 1 div ( D v Vf ) - 2VA” 1 div (l> 2 fL^PD v V /) 

+ VA" 1 div {DlVf) - VA" 1 div (. D 2 fL^PD 2 vv Vf ) (5 . 5) 

+ VA" 1 div ( D 2 fQ(V v v )) - VA" 1 div ( D 2 fL^PD 2 fQ(V v v )) 

+ VA; 1 div (Wo ((W) _1 - id)) - VA; 1 div ( D 2 fL~ x x P (VUiDy)- 1 - id))) 

= - Wo - VA" 1 div (\/ v v). 

In the above, the terms in VA” 1 div appear upon writing Q explicitly. The operator A” 1 o div is 
given by 

A" 1 div (w) = g , 

where g solves 

{ A g = div(io), in fl, 

^L = (w,v), ondn. 

Notice that A” 1 o div is defined up to an additive constant, so VA” 1 o div is defined uniquely. 

Remark 5.2. We notice for further reference, that in (5.5) the first term in every line, except for 
the last and the next-to-the last lines, is linear in /, with the remaining terms being non-linear (in 
/)■ 
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5.2. Local coordinates. In order to have a more explicit description of the operator F acting on 
/, we employ local coordinates. 

Working locally, we choose coordinates (x 1 , x 2 , x 3 ) near <912 such that the domain and its boundary 
are given by 

12 = {x 3 > 0}, <912 = {x 3 = 0}, 


so that 


are tangent to <9f l. We write x 
by the matrix 


<9 

dx 1 

(x', x 3 


and 


a; 3 =0 


d 

dx 2 


£ 3 =0 


In these coordinates, the Euclidean metric is represented 


9 = {gap), a, (3 = 1,2,3, 

with the induced metric on <912 being simply 

9 ij(x', 0 ), i,j = 1 , 2 . 

Also we can assume that <733 = 1 and < 7*3 = 0, i = 1,2. 


Notation 5.3. Unless stated otherwise, Greek indices will run over 1,2,3 and Latin indices over 
1,2. (a/3) means symmetrization on a, /3, i.e., = |(f a/ g + tp a )- The summation convention is 

assumed throughout. 

If U is the coordinate chart, we always assume x G V CC U so that r/(x) € U and rj(x) € U; this 
is always possible when rj is near the identity, which will be the case of interest below. Let 

r = V\dn > 

i.e. 

r(x') = rj(x', 0), 

and write r = (r 1 ,r 2 ,r 3 ). Analogously we have r. With /3 = (/3 1 , (3 2 , /3 3 ), it follows that 

r a ( x ') = /3“(x', 0) + V/ a o /3(x', 0), 

where 


V/“ = g^dj. 


Notice that since /3(<912) = 912, 


/3 3 (x',0) = 0. 


A basis {Ai, A 2 } for the tangent space of r(<912) = 07/(12) is given by 


Xi = dir = Dr ( 

ox 1 


£ 3 =0 


where Dr is the derivative of r. Component-wise, 

= g^d^f o pdi p v + d^d^f o p + sjdip, 
where we have used /3 3 (x', 0) = 0 and 6 ^ is the Kronecker delta. 


The unit (inward) normal to r(<9!2) is 

N a = 


,-a yP y'l 

£ d'yX 1 A 2 


^/e\ u e x<TT X?X°XZXT 
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where e a p^ is the totally anti-symmetric tensor (with the convention £123 = 1). The metric g induced 
on r(d£i) is 

g ij = g(X i ,X j )=g af} X?X?. 

More explicitly, 

9 ij = o pd aX f o pdiP v d 3 p x + 2 d kv $ o pd {i p v d 3 ) p k 

+ 2 d av f o pd^f o pd {i g aiJ d 3) p v + 2g ak d lx f o pd^ i g afl d j) p k 
+ 9ap di9 a>x djg^ x d^f o pd x f o p + g k idip k djP l , 

where (ij) means symmetrization in i,j (see notation 5.3). The second fundamental form of drj(£l) 
is defined by the equivalent expressions 

Aj = —g(S7iN, Xj) = g(N,ViXj), 

where V is the Levi-Civita connection, and the negative sign on the first equality occurs because N 
is the inner normal. Component-wise 

A, = g^N^tXf 

= SafSiX^N" + g a Ai x f N “. 

where T 7 ^ are the Christoffel symbols. Computing get 

diXf = g^d^f o pd t p v d 3 p x + g^d^f o pdijP" + 2 d {i p v d j)g ^d^f o p 
+ d l 3 g ail d t J o p + 5 a k dijP k . 

Also the mean curvature is dehned by 

A = g ij Aj. (5.7) 

Note that the pressure splits into an interior term and a boundary term, p = po + kAh (see (4.3)), 
and since (4.12a) involves XAh 0 V rather than V Ah 0 V, we shall need expressions for quantities on 
the boundary r(d£l) = drj(Q). This amounts to setting P = id in the above expressions, leading to 

X? = g°' J d ia f + d^rd.j + <5“ , (5.8) 


diXj = g^d^ijf + 2d (ig^dj^f + d^g^d^f. 


(5.9) 


~ Q _ g a3 + T“(f) 

x /l + 2g^T x {f) + Tx{f)T\f) 1 

and 

9ij = g^d^fdvjf + 2dijf + 2d lx{ ifd j) g tlu d v f + g a .pd i g (X,l d lL fd j g fiv d v f 
+ gajdig^dpf + gpidjgP^d^f + g i3 , 


where g is the induced metric on r(d fi) and 

T a (f) = eafr^d^f + d 1 ^d li f)(g^d ia f + d 2 g^d ll f) 

+ ta^ig^d^if + dig^d^f) + eai 7 (3 7M 9 At2 / + d 2 g’ ytl d l _ l f). 


(5.10) 


(5.11) 


(5.12) 


Recall that 533 = 1 and g^i = 0. Had this not been the case, the first term inside the square root in 
the expression for N would be g 33 . 
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We shall also need the inverse of the induced metric g, which is given by 


r' = 


i 


522 “512 


, _ . - - . • (5.13) 

detg V “521 <7n 

We proceed to obtain more manageable expressions than those above. Using Taylor’s theorem 
with integral remainder, 

1 




= 1 - ( 5 A 3 7a(/) + -T A (/)T A (/) 


1 -t 


[1 + t(2g X 3 T\(f) + Tx(/)T A (/))] 5 


■ dt 


(5.14) 


= 1 + M(f), 

where A4(/) is defined by the above expression. Combining (5.10), (5.12), and (5.14), 

N a = (g a3 + T a (f))(l+M(f )) 

= g a 3 +M a (f), 

where M“(/) is defined by the above expression. 

We write (5.11) as 

9ij = 9ij T A/)j(/), 

with Mij(f) defined in an obvious way. From this, 

det( 5 ) = det ( 3 ) + £>(/), 

where 'D(f) contains all terms in det(g) that depend on /. Using Taylor’s theorem again, 


(5.15) 


(5.16) 


1 


1 


_ = _ g(/) 

det ( 5 ) det {g) (det (g )) 2 


1 - t 


1 + t 


v(f) 

det(9) 


dt. 


(5.17) 


From (5.13), (5.16), and (5.17), we find 

/ „ 


r l = 


_ V{f) 

det ( 5 ) (det(sO ) 2 


1 -t 


1 + t 


P(/) 

det(9) 


■ dt 


522 + A/2 2 (/) —512 - A/l2(/) 

-521-A/21 (/) 511+A/’ll (/) 


= 9~ 1 +Hf), 


where the matrix J-(f) is defined by this expression and contains all contributions in /. We write 
the above component-wise as 


5 " 5 iJ • r-'(j). 


(5.18) 


We remark that — g a 3 d a f = d v f (the outer normal derivative of /). Also, since the boundary 
Laplacian A is given by (see remark 5.3) 


A = 


1=9. (VWw^i ), 


V\g\ 
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with |< 7 1 = det(gij), we can write 


g ij 9 a3 d ija f = -Ad u f + X(f), 

where X(f), which is defined by this expression, involves at most second derivatives of /. We also 
note that in the present system of coordinates, the second fundamental form of dfi is simply T 3 -, 
where T are the Christoffel symbols, and the mean curvature of <9fi, which we denote AgQ, is g %3 T 3 -. 
From these observations, (5.6), (5.7), (5.8), (5.9), (5.15), and (5.18), we conclude that 

F = A o rj = -A d v f - - 2 A m Af + Q (3) ( d 2 f , d 3 )f + QW ( d 2 f , d 2 )f + A d n on dSl, (5.19) 

where Q® (3 2 f, d 3 ) and Q ty2 \d 2 f,d 2 ) are, respectively, third- and second-order pseudo-differential 
operators. In the present coordinate system, they take the following form 

Q (3) (<9 2 /, d 3 )h = {r j (f)g a3 + g ij M a (f) + r\f)M a (f)) d aij h, (5.20) 

and 

Q {2) (3 2 f, d 2 )h = (F 3 (f) + g 13 M a g a fj + g af) F j (f)M a (f)) {2d {l g i >‘d j) h + d^d^h) 

+ (r %r\f) + g iJ g«AlM a (f) + g aP rf^(f)M a (f)) (g» x d Xj h (5 . 2 i) 

+ 3A X d x h) + Q ij (h). 

In the above, Q 33 (h) is the linear operator in h with /-dependent coefficients, naturally associated 
with the term g a AijX 13 + T 3 jT 13 [f) + g a gg l3 T 3 jM a (f) that figures in the mean curvature. 

More precisely, from our definitions it follows that 

g a A^r 3 (f)M a (f) + T%r\f) + gapg ij r 3 jM a (f) = a\f)d i3 f + ^'(/)%/ + c a (f)d a f , 

where a 1 , b 13 , and c a are smooth functions of / and its derivatives of order at most two, provided 
that / is small (see the beginning of this section). Then, 

Q ij (h ) = a\f)3 a h + b i3 (f)dijh + c a (f)3 a h. 

The explicit form of of Q 13 is too long and cumbersome, and will not be necessary for our purposes. 

Summing up, (d 2 f, d 3 ) and QA{d 2 f,d 2 ) are, respectively, third- and second-order pseudo¬ 

differential operators whose coefficients depend smoothly on / and its derivatives of at most second 
order, and such that Q^ 3 \d 2 f,d 3 ) = QA(d 2 f,d 2 ) = 0 if / = 0 (in particular, Q^ 3 \d 2 f,d 3 ) and 
Q( 2 \d 2 f,d 2 ) contain no zeroth order terms in /) 


6 . Analysis of f\ dn 

We shall now work modulo constants. This suffices to our purposes since we are interested in 
obtaining estimates for V/. Doing so, we can drop the gradient in front of every term (5.5), obtaining 
an equation for / which, upon restriction to the boundary, gives an equation for /It reads, after 
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using (5.19), 


/- 

+ 

+ 


KAd u f - i 

div | 
Q { 2 ) (d 2 f,d 


KAanAf + kQ® ( d 2 f , a 3 )/ + kQ® ( d 2 f , <9 2 )/ 
(vW^(-Aft,/ - ^AonA/ + Q { 3 \d 2 f,d 3 )f 
2 )/))(-id +( J £AT 1 )' 


- kA" 1 div [ J D 2 /Lr 1 p((v^(-AA/ - ^A?nA/ + Q (3) (<9 2 /, 9 3 )/ 

+ Q { 2 \d 2 f,d 2 )f)](D?j)-^ 


+ 2 A ” 1 div ( A,V/) - 2 A ; 1 div ( D 2 fL^PD v Vf ) 

+ A ; 1 div (H 2 „V/) - A ; 1 div ( D 2 fL^PD 2 vv Vf) 

+ A ; 1 div (D 2 fQ(V v v)) - A ; 1 div ( D 2 fL^PD 2 fQ(V v v )) 

+ A ; 1 div (V< 7 o ((^)~ 1 - id)) - A ” 1 div (D 2 fL^P (V ®,^ -1 - id))) 
= — A ” 1 div (V v v) , on <9fi. 


( 6 . 1 ) 


Above, H is the harmonic extension operator in the domain dr/(Q) = 5((id+V/)(J))), and we recall 
that T-L^j is given by (see notation 2.3) 

W-rfiti) = (7i(h o rj~ 1 )) o rj, 

for h : dfl — > M. In (6.1), the function h in the argument of is 

-A d v f - \A d nAf + Q ( 3 ) (d 2 /,d 3 )/ + Q^(d 2 f,d 2 )f. 

The operators Q^ 3 \d 2 f, d 3 ) and Q,^ 2 \d 2 f^d 2 ) were defined in section 5.2, although their precise 
form will not be important here. Rather, it will be important that 

|| Q.( 3 \d 2 f, d 3 )h \\ S)d < C || / || s + 2 ,a|| h || s+3)9 , 

and 

|| Q {2 \d 2 f, d 2 )h || Si a< C || / || a+2l s|| h || s+2 ,a • 

Remark 6 . 1 . In ( 6 . 1 ), the terms in nAaci have been dropped because Aao. is assumed constant, and 
therefore these terms do not contribute in (5.5). Had Aan not been constant, such terms, linear in 
a, would grow without bound in the limit k —>• oo, agreeing with the intuitive idea discussed in the 
introduction that no convergence should be obtained in this case. The term —% corresponding to 
—Vgo does not figure in equation ( 6 . 1 ) either in that go vanishes on <912, since po vanishes on < 97 /( 12 ). 

( 6 . 1 ) was derived under the assumption that / is a sufficiently regular and small solution to 
(4.12). In particular, / was thought of as defined over 12. When viewing (6.1) as an equation for 
f\ d £i, it is important to realize that it depends on the way f\ gn is extended to 12. We are ultimately 
interested in the case when such an extension is carried out using proposition 3.1. However, in order 
to apply techniques of continuous semi-groups, we need the extension to be defined for any function 
in .£P + 2 (<9I2), s > | + 2, and not only for those that are sufficiently small. 
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Let if : M —>■ [0,1] be a smooth monotone function such that if = 1 on (—oo, and if = 0 on 


[2y,oo), where Jo is given by propostion 3.1. Given h € // S+ 2 (<9H), consider the problem 



h \\ 2 s+2,dW(f) = 0 


in fl, 
on dQ.. 


( 6 . 2 a) 

( 6 . 2 b) 


where M is as in (3.2). 

To solve (6.2), we argue as in the proof of proposition 3.1. We define a map 


H s+ 2 (dQ) x H s+ i{Q), 


h \\ 2 8+ 2,aWU)), 


F : H s+2 (dn ) x H S+ 2 (Q ) 
by F(h, f) = (/ — h, Af + ■ 

and notice that F is C 1 , since in a Hilbert space the square norm function is continuously differ¬ 
entiable. Moreover, its linearization is given by (3.4). Thus, arguing as in proposition 3.1, we find 
a small r > 0 and a map ip : B s r +2 (dPl) —>• H s+ 2 (H), where ip(h) = f gives a unique solution to 
(6.2). Noticing that (3.4) does not involve if, if Jo > 0 in proposition 3.1 is sufficiently small, we can 

take r > Jo- Solutions to ( 6 . 2 ) agree with those of (3.3) if || h || s+ 2 ) a< and with the harmonic 

extension of h if || h || s + 2 ,a> Defining (p as the harmonic extension when || h || s + 2 ,a> 

we obtain a map ip : H s+2 (dQ) — > Lf s+ 2 (fl), given by <p{h) = /, where / solves (6.2), and agrees 
with the solution of (3.3) if h is sufficiently small. 

Furthermore, since ip is continuous, given e > 0 we can choose Jo so small that || ip(h) || s+ 5< e. 
But if / is a solution of (6.2) with || / || s+ s< e, and e is sufficiently small, then by elliptic theory the 
solution obeys the estimate 


/ || s+ s< C || h ||s+2,9, 


(6.3) 


where the constant C depends only on e, s, H. Finally, if hu is the harmonic extension of h, from 
standard elliptic theory we have the estimate 


f-hn \\ s+ i<C\\f \\ 2 <C\\f\\ 2 +2>d , 

£ 2 


where the last inequality follows by (6.3). In particular 

II f~hH || s+ 


s< C5l 


(6.4) 


(6.5) 


Definition 6.2. We shall call the solution / of (6.2) constructed above the i^-harmonic extension 
of h. 


Notice that the ^-harmonic extension depends on the choice of which, in turns, depends on Jo- 
We fix these quantities once and for all. 

Notation 6.3. Recall that H o(<9fl) denotes the Sobolev space modulo constants, and that Aon is 
the mean curvature of 

Lemma 6.4. //Jo is sufficiently small, where Jo is defined as above, then the operator (which depends 
on So) — A<9j, — bMdoA : H^ + 2 (dfi) C ff(j(dfl) -A- Hq( 8 Q) C Hfjldfi), where s>^ + 2 = | + 2 and 
d v is computed using the fj-harmonic extension to is an elliptic, positive, invertible, third-order 
pseudo-differential operator. 

Proof. If d v is computed using the harmonic extension, then by known properties of the Neumann 
operator [65], —Ad u is an elliptic, invertible, third-order pseudo-differential operator. Thus, the same 
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holds true when the ^-harmonic extension is used because of (6.5), and in particular —Ad u — ^Agn A 
is also a third-order pseudo-differential operator. 

For what follows, recall that the constancy of Agn implies that dQ is a sphere [2], and therefore 
Agn > 0 . _ _ _ 

Suppose that —A d u f — ^AAgcif = 0. Since the kernel of —A is zero in //q + 2 (30), we have 

d v f + ijAgnf = 0. 

Using integration by parts, 

0 = [ f(dJ + Uanf)= [ |V/| 2 - / MU / Wl+2,aW(f) + \ [ -W 2 , (6.6) 

Jan z Jn Jn z Jgn 


where we recall that Af is given by (3.2). Using the Cauchy-Schwarz inequality, (2.2), (3.2), |^>| < 1, 
the interpolation inequality, and the fact that under our assumptions s + f > 6 , we see that 


/ 

Jn 


Mil / lls+2,9 


W(f) 


< C II / Ikll / || 2 || / Hue'll/ 


,i+l, 


4 
I 5 

16 • 


Also since we are working modulo constants: 

II / ||o< C || V/ 11o • 

Combining ( 6 . 6 ), (6.7), and ( 6 . 8 ), we find 

Agn ii r 112 i n t ii2 


0 > 


/ llo + II / Hi (i- II / Hill / III), 


(6.7) 


( 6 . 8 ) 


and thus we conclude that / = 0 if 5q is very small (recall (6.3)). The invertibility result now follows 
from the Fredholm alternative, and positivity from that of Ad v . □ 


Remark 6.5. In what follows, we will use the fact that A v 1 o div is a bounded linear map between 

R 0 S (U) and #o +1 (^)- 


In order to study (6.1), we first consider a linear equation, with /-dependent coefficients, naturally 
associated with (6.1). Since the map / e+ ( Drj)~ l = (id +D 2 f)~ 1 is not linear, we consider the 
following linear map connected to (Drf)~ l . For / small in 77 S+ 2 (U), the Sobolev embedding theorem 
tells us that ( Drj ) -1 is well-defined, and one can write 

(Drf)- 1 - id = (Dr})- 1 {id-Dr}) = -(DTj^D 2 /. 


This suggests considering the linear map 

B f (h) = -(D?j)- l D 2 h, 

which satisfies the estimate 

II B f (h) || s+ i< C (l+ || / || s+ §) || h || s+ |< C || h || s+ |, 
provided that / is small. 

Notation 6.6. It is convenient to write Bf(D 2 h ) for Bf(h ) to facilitate keeping track of the number 
of derivatives. 
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We are thus led to the following linear equation for h: 

h - kAOJi - ^KA an Ah + kQW ( d 2 f , d 3 )h + kQ® ( d 2 f , d 2 )h 


+ kA,. 1 div 


V7^(-A d u h - ^A dn Ah + Q^ 3 \d 2 f,d 3 )h 


+ Q {2) (d 2 f, d 2 )h )) (- id HD^n)- 1 ) 

- kA~ 1 div \D 2 fL^ l p[(yH^-Ad v h - i A m Ah + Q^ 3 \d 2 f,d 3 )h 
+ Q^(d 2 f,d 2 )h)](D7j) 


\-i 


(6.9) 


+ 2A” 1 div ( D v Vh ) - 2A” 1 div ( D'*fL^PD v S]h 


-l 


l2ir—1 


+ A" 1 div {D 2 v Vh) - A" 1 div (D 2 fL^PD 2 vv Vh) 

+ A~ l div ( D 2 hQ(V v v )) - A” 1 div ( D 2 fL^PD 2 hQ(V v v )) 

+ A” 1 div ( Vq 0 B f (D 2 h )) - A ; 1 div (D 2 fL^ 1 P ( Vq 0 (Bf(D 2 h ))) 

= —A ” 1 div ( V v v ), on <912. 

In (6.9), it is still understood, as before, that rj = id +Df. A simpler, also linear, problem connected 
to ( 1 . 1 ) and ( 6 . 1 ), was studied by the first author in [18]. 

We shall write (6.9) as a first order system. In view of lemma 6.4, the operator L : I7g + 3 (cAI) —>• 
Hq( 8 Q) given by 

£ = -Ad u - -Aan^i 

s-\~— 

(where d u is computed using the ^-harmonic extension to II) has a square root S : H 0 2 (<912) —>• 
Hj&dn), i.e., 

S 2 = £. (6.10) 

(see, e.g., [36]). Lemma 6.4 also implies that exists. 

Letting z = (y/iiSh,h) we will construct solutions to (6.9) by analyzing the following system for 
z: 

dtz + A K (t)z = Q, onfAI, (6-11) 

where 

g e c 0 {[o,T},H s 0 + 1 2 (dn) x h s 0 + ^ {dn)) n cWo^}, h*- 1 ^) x H^-\dn)), 


T > 0, is given by 

With 


Q = (0, —A v 1 div (V v v)) 


( 6 . 12 ) 


/ € L°°([0,T],H s+2 {dn))nC o {[0,TlH s+ \dn)), 

and / extended to 12 via the ^-harmonic extension, 

v € C°([0, T], 2/ s (12, R 3 )) with div(u) = 0, 

and 

q o eC°([0,T],H s+1 (n,] 
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s > | + 2 , A(t) = A K (t) is a one-parameter family of operators A(t) : H 0 2 (dQ) 
as a finite sum 

12 

A = '£ / Z i , 

|*o 

with the operators Z % = Zi(t ), in turn, given as follows. 

z 0 = ^(° 1 —q'\ s. 


Hq 1 ((?n) written 


where 


Zi = V^(\ o')Wi» 


m : H S 0 + 2 ( 8 Q ) -A 
h ha (Q ( $(aPf,& i )S- 1 h 


an 


Z 2 = V ^(J °W 


where 


w 2 : 5 (an) -a // 0 S (5Q), 

h ha (q ( 2) (5 2 /, a 2 )5^'/i) 


= o)w 3 , 


where 


w 3 : fl^ +3 (an) -+ h ( ]-'(5Q), 

h ha (A ; 1 div (v7l^(Sh + Q ( 3 ) (9 2 /, d 3 )S _1 /i + Q ( 2 ) (<9 2 /, S 2 )^" 1 /^ (- id +(£>t/) 


-In 


an 


Z 4 = -v^fj 


where 


W 4 : fl'o +a W -»• AT 1 ^). 

h HA (a ; 1 div + Q { 3 \d 2 f,d 3 )S~ 1 h + Q { 2 \d 2 f,d 2 )S~ l h )) (Df /)" 3 

'0 o' 


an 


where 


2 , = 2 


0 1 


Wo, 


W ,5 : ^o + 5 ( 9 f 2 ) -a fl-* 5 ( 9 f 2 ), 
h ha (A " 1 div (D v Vh)) | 9n . 
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where 


where 


where 


where 


where 


^ 6 =-2 



W 6 , 


W 6 : H* + * (Oil) -A H'l ^(dil), 
h ha (A" 1 div (D 2 fL^PD v Wh))\ dn . 


Z 


= — (° 
^V 1 °y 


Wr, 


W 7 : H^(dQ) -A H(]{dn), 

(A" 1 dr viDlVS^h))^. 


Zs 



m, 


W 8 : H^(r)n) -a tf 0 s (<9C2), 

h* {AZ l d\v{D 2 fL^ 1 PD 2 vv VS- 1 h))\ dil 


Zq — 


0 0 
1 0 


W 9; 


w 9 : ( 90 ) -a h° + * (dn), 

h ha (A; 1 div((Z? 2 5~ 1 h)Q(V„i;)))| ao . 


2 io 



Wio, 


Wio : flS +i (sn) -A flj + »(sn), 

/lHA (A” 1 div ((D 2 fL^[ l PD 2 S~ l h)Q(y v v 


z 


= — (° 
v^V 1 


Wn, 


Wii : -A i?o + 2(c)0), 

/lHA (A” 1 div(Vg 0J B/(-D 2 5" 1 /i)))| 9o . 


^12 



W 12 , 


an • 


where 
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where 

W12 : -A ^(dH), 

h HA (A; 1 div ( y D 2 fL^ 1 P (WUBf(D 2 S- l h))))\ dQ . 

In the above, we use the fact that div(w) = 0 implies Q(S7 v v ) € H s {Ll). In these expression, again 
we assume that quantities are extended to II via the '0-harmonic extension when necessary. 

Notation 6.7. We denote by M. = ■Mr/,t,s +2 the set of functions / : [0, T] -A I7q + 2 (<9II,M) such 
that 

II fit) || s+2 ,c>< R , and || /(f) - /(t') || s+ i )S < ^|f - t'|, 

0 < t,t' < T, where s > ^ + 2. 

Remark 6.8. For future reference, we note that / € C°([0,T], H s+1 (dLl)) if / € At. Indeed, an 
application of the interpolation inequality gives 

II fit) ~ fit') || a +l,0 <|| fit) ~ fit') II s 3 + i 9 II fit) - fit') 111+2,9 
< (2^ 2 R)3|t-t / |l. 


Remark 6.9. In the proofs below, we make extensive use of (2.2) to estimate the several products 
involved. We do not mention the application of (2.2) at every step in order to avoid repetition. 

Notation 6.10. In the ensuing estimates, it will be important to keep track of the dependence of 
several constants on the constants Kq, K' 0 and that appear in the hypotheses of the statements. 
This is done by writing C = C(Kq) and similar expressions. The dependence on fixed quantities 
such as s, II etc, however, will not be indicated, nor will be the dependence on %, q 0 , v, and v, which 
are fixed throughout the theorems of this section. 


Proposition 6.11. Assume that 

veC° (|0,T],R s (H,M 3 )) , div(w) = 0, 

and 


q 0 €C°{[0,T\,H a+1 (Sl, 


s > | + 2, T > 0. Let f G M. satisfy 


f II < Ko 
J |U+2,<9< 

V K 


for some constant Kq, 0 < t < T, and assume that f is extended to II via its if-harmonic extension. 
Finally, let 

G € C 0 ([0,T],R 0 s+ ^(dH) x Hq + ^ (<9H)) n C' 1 ([0, T], Hf~ l (dLt) x H s Q ~ l {dtl)), 

s+l S+ 1 

and zo € H 0 2 (<9H,R) x H 0 2 (<9H, M) be given. Then, if n is sufficiently large, there exists a unique 
solution 

z € C°([0,T],R 0 s+ ^(5H,R))nC 1 ([0,T],R 0 s - 1 (OII,M)) 

of (6.11), satisfying z(0) = zq. (we note that the operators in equation (6.11) are constructed with 
the help of the ip-harmonic extension.) 
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Proof. Since / is defined on Q by its i/’-harmonic extension, we find that (using (6.3)) 

II / || r+ i< C || / || r ,9, 

and thus 

n m ilj< A 

We consider (6.11) as an abstract evolution equation in X = Hg(dfl) x for the unknown 

2 . Our goal is to verify that the operator —A(t) satisfies the conditions of theorem 2.5, and then 

s+ i s+ i 

to apply theorem 2.6. For this, we take Y in that theorem as Y = H 0 2 (<9D) x H 0 2 (dfl), and let 
Hq (dQ) x Hq (cKl) be the domain of Aft). Denote 

Mt) = A(t) = A= -L A K (t). 

V k 

As / is small in H s+ 2 (H) if k is large, the Sobolev embedding theorem guarantees that ( Drj ) _1 = 
(id +D 2 /)- 1 and Lf 1 are well-defined. Using Taylor’s theorem with integral remainder to estimate 
(id+D 2 /)^ 1 , similarly to what was done in the calculations of section 5.2, we obtain 

|| (DTj)- 1 || r _!< <7(1+ || / || r+l)j (6.13) 

1 < r < s. In the above, and in what follows, we assume that k is sufficiently large to allow us to 
bound powers of norms of / by terms linear (in the norms of) /. A similar application of Taylor’s 
theorem also yields 

|| - id +(Drjy 1 || r _i< C || / || r+ i, (6.14) 

1 < r < s. From (5.20), (5.21), and lemma 6.4, we obtain the estimates 

II Wih ||r—1,9 = || Q^\d 2 f, d 3 )S~ 1 h ||r-i,9< c II / || r+ ||| h || r+ i ja , (6.15) 

and 

II mh Hr—1,9=11 Q (2) (3 2 /, d^S^h ||r-i,9< c II / Ilr+|II h || r _i >a , (6.16) 

1 < r < s. 

We claim that if g € I7 s_ 1 (9f2), then PLjj(g) G H s ~ 2 (D). To see this, notice that G = 'Hjj(g) means 
G = G orj, where G solves 

f AG = 0, in 

1 (7 = 5 o jf -1 , on drj(Q). 

But using notation 2.3, this can be rewritten as 

f A?jG = 0 , in S 2 , 

|(7 = g, on dD. 

If G is defined over 0, we compute, 

A^G = {(Dr})- 1 )%(d ce jG)((Drj)- 1 )fr + dadDrfi-^dpG, 
where ((Drj) -1 )^ are the entries of (. Drj ) _1 . Thus, has the form 

A ri G = a ap d a pG + b a d a G. 
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We see that that a Q/3 is in H s+ 5 (fl, M 32 ) and b a is in H s ~5 (J7, M 3 ) since V/ € H s+ § (0). Furthermore, 
as V/ is small in H s+ 2 (f2) and s > | + 2, we know that ( a a ' 3 ) is positive definite. Thus, A^ is an 
elliptic operator, and the solution G will be in H s ~z{£l) if g is in f7 s_1 (<9f2). Moreover, we have the 
estimate 

II 'Hjjig) || s _i< C(f) || g IU- 1 , 0 , 
where the constant C(f) depends on || V/ || s . 3 . 

Combining the above properties of T~Ljj with (6.14), (6.15), and (6.16), it follows that 

II W 3/1 ||)—1,3 


= 11 AF 1 div 


VUrXSh + ( d 2 f , d 3 )S^h + Q (2 ) ( d 2 f, d 2 )S^h)) (- id +(Drj) 


Ir—1,0 (6.17) 


3 || h 
' 2 


<C\\f 

1 < r < s. Similarly, (6.13), (6.15), and (6.16) give 

|| A” 1 div ^D 2 fL^p{vHrj{Sh + Q {3) (d 2 f, d^S-'h + Q {2 \d 2 f , d 2 )^ 1 /*)) (. D7j)~ l 

<C\\ f || r+ 3|| h || r+ i ia> 
so that 


r—l,d 


mh || r _i,a< C II / 


1 r+i 


h 


lr+ia> 


(6.18) 


if 1 < r < s. 

If k sufficiently large, (6.13), (6.15), (6.16), (6.17), (6.18), and the form of A, imply that A is 
arbitrarily close to 

Zq = —j= z $ 

V K 

in 7?(i7g(9fl) x Hq(80,),X), with | < r < s + thus in particular in B(Y,X ) and in B(H 2 (dfl) x 
Hz(dfi), X). Indeed, (6.15), (6.16), (6.17), and (6.18) show that 




is close to Zq, and the remaining terms in A are lower order operators multiplied by 

Now we are ready to verify the several hypotheses of theorem 2.5. We start by first showing that 
theorem 2.4 can be applied. By the operator A, we in fact mean its maximal operator, so that A is 
in fact closed in X. We shall not, however, make notational distinctions between A and its extension 
(see [33] for details). 

We start by showing that A generates a semi-group. Notice that Zq + A, A € R, is invertible (with 
bounded inverse), and so is A + A. Write 


A — (A — Zq) + Zq — £ + Zq. 

From the above estimates, we know that 

II Sz ||o,a< C(K 0 )s || z ||| i(9 , 

where e can be as small as we want provided that k is sufficiently large. Set 

w = {A + A)z, 


(6.19) 

( 6 . 20 ) 


( 6 . 21 ) 
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and use the fact that Zq is skew-symmetric to obtain 

O, z ) 0 ,a = (£z, z) 0 ,d + |A| || z ||o,9, 

from which we get (using the Cauchy-Schwarz inequality and (6.20)) 


|A| 


\d <11 2 11 0,0 1 1 w ||o,a + II £ z Ho, dll ^ 11 o,c? 


<11 2 ||o,9ll W ||o,9 +C(Kq)s II z ||3 ia || z ||o,9 • 

Next, we show that 

II ^ ||§, 9 < C(|A| || z ||o,9 + || w || 0 ,a). 

Recalling the definition of Z 0 , (6.19) and (6.21), we can write w more explicitly as 


ici 

W2 


From the hrst row of (6.24), 
which implies 


A -5 
S + £21 A + £22 


Z 2 = A S l zi — S 1 w\ 


( 6 . 22 ) 

(6.23) 

(6.24) 


(6.25) 


II -2 ||§, 9 < +C(|A| || zi ||o,9 + || wi ||o,a). 

The second row of (6.24) gives 

Zi = -S- 1 £ 2 1Z 1 — A S~ 1 Z2 ~ S~ 1 £ 2 2Z2 + <S~ 1 W2, 

which gives (with the help of (6.20) and (6.25)), 

II Z! II3 g< C(K 0 )e II Zi II3 g +C|A| || z 2 || 0 ,c> +C(|A| II Zi ||o,a + II W! Ho ,9 + II W 2 Ho,a)- (6.26) 
2 ’ 2 5 

The term C(Kq)£ || zi |j 3 g can be absorbed on the left-hand side by making e small, and then (6.25) 
and (6.26) combine to give (6.23). 

From (6.22) and (6.23) it now follows that 


IA | || z ||o ,0 < C(Kq){\ + e) || z ||o,d|| w ||o,9 ~\~C{Kq)e\ 


I o,9 


Thus, if e is sufficiently small, we find || z ||o,9< 


C{K 0 ) 

|A| 


w 11 0 , 9 j ° r ? since z = {A + A) 1 u;, 


(-4 + A) w ||o,9< 


C(K 0 ) 

|A| 


w ||o ,a 


(6.27) 


From theorem 2.4, we conclude that the family A = A(t ) generates a continuous semi-group, and 
that (i) of theorem 2.5 is satisfied. 

We now consider (ii) of theorem 2.5. As shown above, A is close to Zq if k is large. Thus, 
A = (A(t))~ 3 "" gives an isomorphism from Y to X, and AA(t)A -1 generates a semi-group on X. 
This implies, by the results of [33], that e rA ^ t ' > restricts to a semi-group on Y. Moreover, setting 

t (zr,z 2 ) = ((A(t))^z 1 ,(A(t))^z 2 ) 0 ,9, 

one obtains an inner product that generates the topology of Y. Then, an argument similar to the 
one leading to (6.27) gives 

C(K 0 ) 


(A + A) ||o P (t)< 


|A| 
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From this, the remaining conditions of (ii) follow. From the continuity of 1 1 —> f(t ) in the H s+l (dPl) 
norm (remark 6 . 8 ), the assumptions on v and % and and the form of A, we obtain the continuity of 
t eA A{t) in B(Y,X). The remaining hypotheses of theorem 2.5 are routinely checked, and therefore 
we obtain the desired evolution operator for A, and hence for A as well. 

Invoking theorem 2.6 and given an initial condition, we obtain a solution to (6.11), namely 

zGC°(lO,T],Y)nC 1 (lO,T],X) 


It remains to verify that t H > A(t)z(t) is continuous with respect to the H s 1 norm. Because 
z € C 0 ([0, T], Y), we have D 3 S~ 1 z\ dii € C°([0, T], iFg - 1 (c>f2) x Hq~ 1 (8Q)). The result now follows 
since the coefficients of A depend on at most two derivatives of /, and are therefore continuous in 
the FP - 1 (<9I2) norm because / G A4 (see again remark 6.8). □ 


Proposition 6.12. In proposition 6.11, assume further that 

v € C° ([0,T\, H s (n,R 3 )) n C\[0,T], H s ~% (n,R 3 )), 

q 0 € C° {[0,T},H s+ \n,R))nC\[0,T],H s -Hn,m, 

and 

f G C'°([0,T],FT 0 s+ 2 (5O,M))nC7 1 ([0,T],Fr 0 s+ ^(5O,M))nC 2 ([0,T],i? 0 s - 1 (aO,M)) 


satisfies 


II / lls+1,9— K 

for some constant Kq. Then the solution z satisfies the estimate 

II z{t) || a+ i Q < Ki || z( 0 ) || a+ i Q + sup II G(t) ||s-i,a 

2 2 0<r<t 

Ki{l + t) ( 

+ sup --=— || drGlr) || s _i,a, 

0 <T<t V K 


(6.28) 


for a constant K\ given by K\ = sup 0 < T <r V(t), with V a continuous function of 
II v(r) IU II v(t) || a _3, || % (r) || s+l > II 9o( T ) \\ s - 1 ,T,K 0 , and K' 0 . 


Proof. It follows from the calculations in the proof of proposition 6.11 that 

1 


i-i 


A K 1 = —j=A~ , 
C(k 0 ) 


and that the norm of A K 1 is bounded by for large k. Hence 




C(K 0 ) 


W s -l,9 


(6.29) 
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Next, invoking Duhamel’s principle, we have 

rt 

ztt) = ’ 


(t) =U(t,0)z(0) + f U(t,r)Q(r) dr 

Jo 

= U(t,0)z(0) + ( U{t,T)A K {T)(A K (r))~ l Q(r)dT 

Jo 

= U(t,0)z(0) + [ d r U(t,T)(A K (T))~ 1 G(T) dr 

Jo 

= U(t, 0)z(0) + (A^t))-^) - U(t, 0)(^ ft (0))- 1 g(0) 

+ f U(t,T)(A K {T))~ l {d T A K {r))(A K {T))~ l G{T)dT - f (A k {t))~ 1 d r Q{ t) dr, 

Jo Jo 

where U is the evolution operator for — A K , satisfying 

d r U{t,T ) =U(t,T)A K (T) 


(6.30) 


and U{t,t ) = I (identity operator), and we have integrated by parts and used 

d T {A K [r))- 1 = -(A^r^drA^iA^r))- 1 . 

Differentiating A K (t), and using the hypotheses on /, %, and v, one checks, with the help of (6.29), 
that 

|| (. A K (T))~ 1 d T A K (T)(A K (T))~ 1 w || s+ i d < C ^ A °_ A o) || w || s _ i a . (6.31) 

2 ’ V K 

We remark that while the expression for d T A K (r) is rather cumbersome, and thus will not be given 
here, the estimate (6.31) follows essentially by counting derivatives, with some of the less-regular 
terms handled via standard arguments. For instance, v G H s ~2 (O, R 3 ). But v is divergence-free and 
tangent to the boundary, thus div(V 1 ) u) G H s ~2 (fi) and (is,V v ii) = —(S7 v v,v) € H s ~ 2 (dLl), so that 
Q(V v v) G H s ~ §(D,M 3 ). Thus (A " 1 div {(D 2 S~ 1 h)Q(W v v))) € H s ~h(Ll), or 

(A ; 1 div ((D 2 S~ 1 h)Q(V v v)))\ dn g H s ~\dn), 

This is needed in order that (6.31) make sense (see Wg in the definition of A K ). 

Then combining (6.29) and (6.31) with (6.30) gives 

II z(t ) || a+ i d < C{K 0 ,K'q)( II z( 0) || S+ 1 g + sup (1 t- 1 ' II S(t) 11,-i.e 

2 \ 2 0 <T<t V N 

(1 + *) ii n n < \ ll A 

+ sup — -j=- || d T G{T) IU— 1,0 . 
o <T<t V K J 

Note that the constant C(Kq, K() appearing in the above estimate has the desired form (see notation 
6.10). This finishes the proof. □ 

Proposition 6.13. Let hi € (dLl,M), and assume the hypotheses of proposition 6.11. Then, if 
k is sufficiently large, there exists a unique solution 

h G C°([ 0 , T], H s 0 +2 (dn, R)) n C^flO, T], (dQ, R)) n C 2 ([ 0 , T], H*- 1 (dn, R)) 


of (6.9), satisfying h{ 0) = 0, and h{ 0) = h\, where where h is extended to Ll via its if-harmonic 
extension. If in addition the hypotheses of proposition 6.12 hold, then there exists a constant > 0, 
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such that, 


, , N „ K 2 „ ; „ AVl+T) 

h(t) lls+2,9 < II M0) || s _|_l ! 0 H — : 

K 2 (l + T) 


II ^ lls+i,9 — II M°) IL+i,a + 
and || h || s _i j( 9 < K 2 , 0 <t<T, where K 2 = sup 0 < r <T 'P(r), w ith V a continuous function of 
II v(t) || s , II v(t) H s _|, II q 0 (r) || s+ i, || %(t) \\ s _i,T,K q , and K' 0 . 

Proof. The argument is more or less standard, so we shall go over it briefly. Let G be given by (6.12). 
Consider the solution to (6.11) with initial condition z(0) = (0,h(0)). Define h as the solution of 

d t h = z 2 , 

with initial condition h{ 0) = 0. Then 

\fHSh(t) = \[k [ d t Sh(r) dr 


= \f~K / Sz 2 (t) dr 

Jo 

= 31(t), 

after using the first line of (6.11), i.e., dtz\ = \[kSz 2 . By inspection we see that h = solves 

(6.9) and has the correct regularity. 

From the properties of S that follow from lemma 6.4, 

(6.32) 


II h ||s+ 2 ,a<|| y/USh || s+ i< Cy/n || h || s+ 2 ,a, 


which combined with (6.28) and the definition of z gives 

1 


h |U +2 ,a < C{K X ) 


• . „ (1 + 1 ) 
.. MO) ll s +i,a + sup 

K 2 0 <r<t K 


G(t) || s _i,a 


(! + *) II a nr \ n A 
+ sup - || O t G{t) |U_ 1 ,c> , 

n<rT^t K / 


0 <r<t K 


and 


h || s+ i ia < C(/C 1 )( || MO) \\ s+ i td + su^^S 


0 <T<t 


G(j) |U—i,s 


+ sup 

0<r<t 


(1 + t) 


K 


d T G(r) || s _i )9 


where K\ is the constant appearing in the conclusion of proposition 6.12, and we have used h( 0) = 0. 
The estimate for h now follows by solving for h in equation (6.9) and using the above estimates for h 
and h. We then invoke the definition of G to obtain the result. Here, we need to check that G indeed 
has the correct regularity. This follows from the fact that v and v are divergence free and tangent to 
the boundary, so that Q(X7 v v ) € H s (Gl,M. 3 ) and Q(V v v) € H s ~ 2 (D,M 3 ) (see similar discussion right 
after (6.31)). □ 
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The main result of this section, theorem 6.16 below, uses the Schauder fixed point theorem and 
the Arzela-Ascoli theorem, which we state for the reader’s convenience. 

Theorem 6.14. (Schauder fixed point) Let S be a closed convex set in a Banach space, and let T be 
a continuous mapping of S into itself such that the image T(S ) is pre-compact. Then T has a fixed 
point. 

Proof. [48] or [27]. □ 

Theorem 6.15. (Arzela-Ascoli) Let X be a topological space and ( Y,d ) a metric space. Consider 
C°(X,Y) with the topology of compact convergence. Let F be a subset of C°(X,Y). If F is equi- 
continuous under d and the set 

T x = {/(*) |/eJ} 

has compact closure for each x € X, then IF is contained in a compact subspace of C°(X,Y). 

Proof. [45], □ 


Theorem 6.16. Assume that 

v€C° ([0,T},H s (n,R 3 )) nC ,1 ([0,r],,H' s -i(fi,R 3 )), div(u) = 0, 


and 

q 0 € C° ([0,T},H s +\n,R)) n C\[0,T\, (Q,R)), 


> § + 2, T > 0. 


Let fi € H 0 


(<9f l) satisfy 


II /. IUj,»< 7 =. (6-33) 

for some constant K 3 . Finally, let Q be given by (6.12). Then, if n is sufficiently large, there exists 
a T' € (0, T] and a solution 

fec\[o,TlH s +\dn,R))nc\[o,TlH s 0 + ^(dn,R))nc\[o,T'},H s 0 -\dn,R)), 

satisfying (6.1) with initial conditions /(0) = 0, and /(0) = fi- In (6.1), it is understood that f is 
extended if-harmonically to Lt. Furthermore, f obeys the estimate 


f 


s+2,d< 



f ILii =< 


K 4 


s + 2 ’ d ~ WK' 


and 

II / ||s-i,a< AT 4 . 

The constant K 4 is given by K 4 = sup 0<r<T V(t), with V a continuous function of 
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Proof. Consider the set Al from notation 6.7 with the metric 

d(f,g)= sup || f(t) - g(t) ||o, d • 

0 <t<T 

We start by noticing that At is a complete metric space in the metric d. Indeed, if / n —> / in the 
metric d, then, from the facts that Br( 0) C H s+2 (dGl) is weakly compact where -Br(O) denotes the 
ball about zero of radius R in the metric d and that the embedding H s+2 (dGl) C H°(dGl) is compact, 
we conclude || / || s + 2 ,a< R■ By the same reasoning we also find that, for each t, f(t ) G H s+ z{d fl), 
and f n {t) —> f(t ) in H s+ ? (dfl) since H s+2 (dfl ) C H s+ z(d£l) compactly. But, 

II f{t) ~ f(t') || s+ l >a <|| f(t) - f n (t) || s+ l )9 + || f n (t) ~ fn{t') || s+ l >a + || fn(t') ~ f(t') || s+ l >9 
<11 fit) - fn{t) || s+ i )S +£\t - t'\+ || f n (t') - f(t') || S+ I . 

Passing to the limit we conclude that || f{t) — f(t') || s+ i a < i\t — t' |, thus f € Al. Therefore, Al is a 

closed subset of the Banach space L°°([0, T], H® (<9fl)). One immediately checks that At is convex. 
Given / € Al satisfying 

K 0 


f l|s+2,9< 




(6.34) 


proposition 6.13 yields a solution h to (6.9) if k is sufficiently large. We shall show that, if T is 
small, k large, and R and £ are suitably chosen, the association / ha h defines a map <!> : Al —^ 
L°°([0, T], Hq (50)), which (i) takes At into itself, (ii) is continuous, and (iii) has image 4>(A1) pre¬ 
compact in L°°([0,T], Hq( 8 Q )). 

To fix the constant Kq in (6.34), we put Kq = K 3 . Choose R = R(Kq,k) = so that (6.34) 
holds for / e At. Thus, in light of proposition 6.13, we obtain a map 


s+i 


C"([0,T],H^ 2 (dn)) n c L ([o,T\,H 0 2 (dV)) n C 2 ([0,T},H°-\dQ)), 

given by <!>(/) = h, where h is the solution of (6.9) with initial conditions h( 0) = 0, h( 0) = f±. Notice 
that <h can be viewed as a map from Al to L°°([0, T], Hq (<9fl)). 

Letting z = (y/nSh, h ), where S is given by (6.10), z satisfies (6.11), and is given by 

z(t) = U{t, 0)z(0) + f U{t,T)G{T)dT, (6.35) 

Jo 

where z( 0 ) = ( 0 ,/i), U is the evolution operator associated with — A K (t ) (see proposition 6 . 11 ), and 
G is given by ( 6 . 12 ). From this and (6.32), 


h 


s+2,9 < 


C(K 0 ) 

\[k 


K 


C(K 0 )T 

fi || a+ i a +- y=~ sup || G(t) || a+ i a 

2 V K 0 <r<T 2 


£( r ) L+i,a» 


z C(K 0 )K 3 , C(K 0 )T 
<- 1 - 7 =— sup 


(6.36) 


0<t<T 


using (6.33). We note that C does depend on I\q because of the estimates for A K (t) in the proof of 
proposition 6.11. Choose T small and k large so that 


C(K 0 )T 


sup 

0 <t<T 


G{t) 


R 

'*+bd< 2’ 


and 


C(Kq)K 3 r 

k < 2 ' 








SURFACE TENSION 


35 


are satisfied. 

Note that this last inequality is possible despite the dependence of R on Kg and k because R is 
of the order of while the left hand side of the last inequality is of order We conclude that 

II h ||s+ 2 ,< 9 < R- (6.35) also gives 


h \\ a+ i 8 <C(K 0 ) || /i || s+i a +C(K 0 )T sup || S(t) 

2 ’ 2 ’ 0<r<T 

< C ( R yl h3 + C(Kq)T sup || g { r ) || S+ 1 g . 

'K o <t<T 2 


h+2,9 


(6.37) 


This implies || h(t) — hit') || s+ i d < £\t — t'\ if we choose £ large enough as to have the right side of 
(6.37) less than l. We conclude that h € M, i.e., <h maps M. into itself. 

Next, we study the continuity of <L. This requires estimating h — h = <J>(/) —<J>(/), where /, / S A4. 

Let z be as above, and ~z = {y/USh, h ), with — A K the operator in (6.11) with / in place of /, and U 
the corresponding evolution operator. We have the estimate: 

II z-z ||o,a < C(K 0 ) [ || (A(t) - A(t))z(t) || 0 ,a dt. (6.38) 

J o 

We recall how (6.38) is obtained. Computing we find 

d s ( U(t, s)U (s, r)y) = d s U(t , s)U(s, r)y + U(t, s)d s U(s, r)y 

= U(t, s)A(s)U(s, r)y — U(t , s)A(s)U{s, r)y. 


Then integrating between r and t, we get 

U(t,r)y-U(t,r)y = f U{t, s)(A(s) - A(s))U(s,r)y ds. 

J r 

Setting r = 0 and y = z(0) in (6.39) gives 

U(t, 0)^(0) — U(t, 0)z(0) = f l4(t,s)(A(s) — A(s))U(s, 0 )z( 0 ) ds, 

Jo 

Setting y = Q{r) in (6.39) and integrating in r yields 

[ {U (t, r) — U(t, r))Q(r) dr = f f U(t, s)(A(s) — A(s))U(s,r)Q(r) ds dr 
Jo Jo Jr 

rt 


( U(t,s)(A(s) — A(s)) f U(s,r)G(r) dr ds. 
Jo Jo 


(6.39) 


(6.40) 


(6.41) 


(6.40) and (6.41) imply (6.38); see [35] for details. 

We need to estimate the difference (A(t) — A(t))z(t ) in (6.38). For this, we point out that the 
operators Zq, Z§, Z 7 , and Zg that figure in the definition of A do not depend on /, being therefore 
the same for A(t) and A(t). Hence, they cancel out in the difference A{t) — A{t). 

Continuing we have 


{A{t) - A(t))z ||o,a = ((A(t) - A(t))z, (A(t) - A(t))z) 0 , d 

= ((A(t) - A(t))* (A(t) - A(t))z, z) 0 ,d (6-42) 

<|| (A(t) - A(t))*(A{t) - A(t))z ||_ r ,a|| z || r ,a 


where {A{t) — A(t))* is the adjoint of A(t) — A(t) in H^dfl). A* and A* are pseudo-differential 
operators or order | whose coefficients depend on at most four derivatives of /, and in the last step 
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we used the generalized Cauchy-Schwarz inequality (a,b)o <|| a ||_ r || b || r , with || • ||_ r denoting the 
negative Sobolev norm, and r < s + b a number that will be conveniently chosen. 

Since (A(t) — A(t))* involves at most four derivatives of /, / G H s+2 (dQ), and s > f + 2 = I ^- + f 
(so that s — 2 > ), we can apply ( 2 . 2 ) to get 


II (A(t) - A(t))*(A(t) - A(t))z ll-^a < C\/~k || <9 4 / || s - 2 ,a|| (A(t) - A(t))z ||_ r+ | 9 

<£W>1|| (A(t)-A(t))z ||_ r+ i 8 , 

VK ' 2’ u 

where <9 4 / symbolically represents terms in at most four derivatives of /, and where we used the fact 
that the s + 2 norm of / gives a - factor. Thus (6.42) reads 


C(K 0 ) 


{A(t) - A(t))z || 0j(9 < — y=- II (A(t) - A{t))z ||_ r+ 3 e || z || S+ 1 Q ■ 

yj Kj 2 2 


(6.43) 


Recall now how (6.11) was obtained from (6.9). Each term in (6.9) is a pseudo-differential operator 
of order at most three, whose coefficients depend on at most second derivatives of /, although this 
dependence may be non-local due to the presence of the operator A ” 1 o div. If p(D)h is one of such 
pseudo-differential operators acting on h, the corresponding term in ( 6 . 11 ) is of the form 

—=p(D) o S~ 1 {y/7iSh) = — =p{D ) o S - 1 zi, 

y/K y/K 


and if p(D)h is one of the operators that acts on h, the corresponding term in (6.11) is of the form 


p(D)h = p(D)z 2 . 


With these considerations in mind, we proceed to the estimates below. Recalling (5.20), we have 

|| V^Q {3) (d 2 ld 3 )S~ 1 z l -V^Q {3) (d 2 f,d 3 )S- 1 z 1 ||_ r+ | |fl 

= yfr II (a ai i(D 2 f,DfJ) - a ai \D 2 f,DfJ))d ai ,S~ 1 z l ||_ r+ 3 >0 
< Cyfc || (a^(D 2 f,DfJ)-a^(D 2 f,Df,f)) ||_ r+ | ia || d aij S~ 1 z 1 || s _i, 9 . 

In the last inequality, we used the fact that z\ € H 0 2 (cKl), so that d a ijS~ 1 zi belongs to -£fQ _ 1 (<9fl); 
we also used the inequality s > | + 2 so that ( 2 . 2 ) is valid. As the coefficients a mj are smooth 
functions of their arguments, we conclude 

|| V^Q {3) (d 2 f,d 3 )S^z 1 ~V^Q {3 \d 2 f,d 3 )S- 1 z l \\_ r+ 3 jd 

< C{K 0 )yfc || f- / ||_ r+ l )S || Z 1 || s+ l j9 

< C(K 0 )y/~K, || / — / ||_ r+ 7 )9 , 

after using 

II Zl lls+i,a— C (6.44) 

since z\ = yfkSh, and || h || s + 2 ,s< ^ 7 =. 

The other terms are similarly estimated so we get 

II (A(t) - A(t))z ||_ r . + | i0 < y/^C(K 0 ) || 7- f \\_ r+ r )d 


(6.45) 
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Then combining (6.43) and (6.45) we get 

|| (1(f) - A(t))z || 0 ,e < C(K 0 )yJ\\ f-f \\_ r+ r i9 . (6.46) 

Here we have used || z || s+ i 9 < C, which suffices since z\ is bounded by (6.44), and || Z 2 || s+ i a 
is bounded by (6.37) and our choice of l. We now choose r = | (which is less than s + 4 since 
s > | + 2), so (6.46) gives 

|| (1(f) - A{t))z || 0 ,a < C(K 0 )y/d(f,f). (6-47) 

On the other hand, invoking (6.32) with s replaced by — 7 and once again using z = (y/nSh, h), we 
get 

|| z — z ||o,a > Cy/~K || h — h ||| )a > Cy/H || h — h 110,0 • (6.48) 

Combining (6.38), (6.47), (6.48), and recalling the definitions of h and h we get 

<*($(/), $(/)) < £VQ T y/d(f,f). (6.49) 

V ^ 

This establishes the continuity of the map <L. 

We now show the pre-compactness of 4>(A4). Recall that / € C°([0, T], Hq (dQ)) if / G M. Let 
{fu} CM,1 an indexing set. Invoking the Lipschitz condition once more, 

II fu(t) - f u {t') ||o,a<|| fv(t) - f v (t') || s+ i )9 < t\t - 1 ' I, 

we see that {f u }v€i is equi-continuous as a family of maps from [0,T] to Hg(dfl). Also, for each 
fixed t, the set 

ifu(t) \u 6 1} 

has compact closure in the i7°-topology in view of the compact embedding Hq +2 (cA 2) C 77g(<9H) and 
the bound || f u || s +2,a^ R- Combined with the continuity of the functions /„, we have thus verified 
the conditions to apply the Arzela-Ascoli theorem, and we conclude that {/„} has compact closure 
in C°([0, T], Hq (<9H)), and therefore in A4 since At is complete. The same then holds for {4>(/ I/ )}^ e x 
because of (6.49), i.e., the continuity of : At —>• At. This shows that 4>(A / f) is pre-compact in 
L°°([0, T], Hq(8 Q)). 

We can now invoke theorem 6.14, namely, the Schauder fixed point theorem, to conclude that T 
has a fixed point in At, i.e., there exists /* € At such that $(/*) = /*. This /* solves (6.1) with 
initial conditions /*(0) = 0 and /*(0) = fi- In view of our choices, we see that /* and /* satisfy 
the further hypotheses on / and / in proposition 6.12 with Kq = K 3 and K r 0 = l. Therefore, by 
proposition 6.13, we obtain the desired estimates for /*, /*, and /*. □ 

6.1. Analysis of / in the interior. The results of section 6 give us a solution to (6.1) with 
appropriate initial conditions. Since the operators in (6.1) involve the ^-harmonic extension (6.2a), 
for / sufficiently small (i.e., k large), the extension of / to H also satisfies (4.12b). We remind the 
reader that (6.1) is (4.12a) without the gradient, or (5.5), restricted to the boundary. 

As before we can drop the gradient in front of every term in (4.12a) and work modulo constants. 
This leads to an evolution equation for / of the form 

if + &t(y,f)(f)+&(v,f)(f)+ < tf(v,q )=0 in Q, 

1/(0) = 0, / = /!, 


(6.50) 
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where srf is a third order pseudo-differential with coefficients depending on v and /, S3 is first order 
with coefficients depending on v and /, and ^ is a lower order operator on v and q. Here, as in 
section 6, q will be a given function which replaces po o rj, and /q and f\ are known functions. 

Let / be given by theorem 6.16. / is then defined on fl and satisfies (4.12b), as stated earlier. 
Since (6.50) gives (6.1) on <9fl, if we plug / into the left hand side of (6.50) we obtain 

f + £?(v,f)(f) + &(v,f)(f)+ e if(v,q) = u, 

where u has the property that is zero in H^(d fi) (recall that we solved (6.1) modulo constants), 
so w is constant on dQ. Therefore, if we work modulo functions that are constant on the boundary 
(which suffices for our purposes), we find that / automatically satisfies the interior equation (6.50), 
and thus (4.12a) as well. 


7. Proof of theorem 1.2: existence 

In this section we prove the existence part of theorem 1.2. (recall remark 1.12). Before doing so, 
we summarize how the argument will be implemented. 


7.1. Overview of the argument. Here we motivate how the iteration yielding a solution to (1.1) 
is implemented, and also fix some notation for future reference, while following the same notation as 
above for quantities that have already been introduced. 

Assume we are given a solution to (1.1), so that rj(t) is a curve of volume-preserving embeddings, 
fj = uo rj, and divti(f) = 0. Write 

u(t) = Pu(t ) + Qu(t ) = Pu(t) + V/i(t), (7-1) 

where h(t) : fi) -A M is harmonic. Recall that u satisfies 

r\ 

-^ + V u u = -Vp, in U (f}xfi(f). (7.2) 

0 <t<T 

Letting Pu = w and using the fact that 

VvhV/i = |v|V/i| 2 , (7.3) 

we obtain, (after applying P to (7.2)), 

+ P(S7 u w + W w Wh) + ViL = 0, in {t}xQ(t), 

0 <t<T 

where we notice that ViL is divergence-free and has normal component equal to (w, (N o rj)'o tR 1 ) + 
(V u w,N). Composing with 77 gives 

, .. ,dw . 

[w o 77) = (— + V U VJ) o 77 

= (— P(V w Vh) + Q(V u w)) o 77 + VH o 77. 

Letting z = w o 77 then gives 

z = Q v ((VuUz)) - P n {{v w ) v (yh o 77)) + vh o 77. 

Next, consider P^(fi). It is a submanifold of H s ( f),R"). Therefore, it has a normal bundle given by 
the L 2 metric on H s (Q,M n ). This metric is of course invariant under right composition by elements 
of D^(fi). A tangent vector to D^(fl) at f3 is of the form v o j3 where divu = 0 and v is parallel to 
<9f2. Hence, a normal vector is of the form V/ o j3. The exponential map from the normal bundle 


(7.4) 

(7.5) 
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to H s (Ll,M. n ) is a diffeomorphism in a neighborhood of Dffi}). Therefore, if 77 is near LF(O), then 
there exist and V g such that 

r/= (id+Vg) o/3. (7.6) 

In other words, decomposition (4.1) holds for all 77 sufficiently close to T> s ^ (11), although it is important 
to stress that in this decomposition, given by the exponential map, g need not to satisfy ( 6 . 1 ) or 
(4.12a)-(4.12b), nor does it need to be in H s+ 5(12) as are the solutions constructed in section 6 . 

Now we can describe the iteration. Assume that we are given a differentiable curve 77 (f) of H s 
embeddings and u G H s (r](Q)). From (7.6) we obtain /5, and thus v, both in H s (Ll). Then we use 
theorem 6.16 to solve (6.1), or equivalently, (4.12a)-(4.12b), for /. Next, we obtain h by solving 

j Ah = 0 , hi 77 (fi). 

S dh _ I lx? i I n vrr I ~-l Tf\ __ 


. ON 


= ((V/ + D v Vf + v) O r] 1 ,N) on drj(£l), 


where rj = id+V/ as usual and N is the normal to 77 (D). Since V/ G H s+ 2 (H), V/ € H s ( 12), 
and v € H S (Q), we find that V/i € H S (Q). Next, solve (7.5), which is an ODE for z. Note that 
V/i G H s (£l) and u € H S (Q) so z € 77 s (D). Finally, let the new 77 be 

77 (f) = id + [ (z + (Vh) o 77 o /3). 

.70 

This also gives the new u by u = z o o rj -1 + V/i. Using our estimates on /, we can then show 
that this iteration has a fixed point. 


7.2. Successive approximations. Here we carry out the fixed point argument sketched above. 
Denote by Emb s (D,M 3 ) = Emb s (D) the space of H s embeddings of 0 into M 3 (not necessarily 
volume-preserving). 


Induction hypothesis (step n, at which the (n + 1 )st quantities will be determined). Assume 
inductively that we are given 

Tin € C 2 ([0, T] , Emb s (0) ) n C 1 ([0, T] , H s (D)) 0 C 2 ( [0, T] , H s ~ § (O) ) , 
satisfying rj n (0) = id and rj n { 0) = uq, where uq is the given initial velocity for (1.1). Suppose that 



II Vn - id || s < R s , 

(7.8a) 


II Vn (0) ||s = || u 0 || s < R 

(7.8b) 


|| Vn Hs— Ri 

(7.8c) 


II Vn{ 0) || s _3< R° s < R, 

(7.8d) 


II Vn IU|< R , 

(7.8e) 

for some constants R s , R), 

and R which will be suitably chosen. Assume also that 



Vn=u n o (id +V/„) 0 /3 n , 

(7.9) 

with u n divergence-free, u n 

€ H s ((id+Vf n ) ofi n (Q)), fi n G i7 s -t((id+V/ n ) ofi n (Q)), 

where fi n is a 


C 1 -curve /3 n : [0,T] —> 'D)(Q) which also satisfies /3 n € H s 2 (D), and 

V/ n GC°([0,T],77 s+ l(D))nC7 1 ([0,T],i7 s (D))nC' 2 ([0,T],77 s -|(D)). 

Let v n be given by f3 n = v n o (3 n . Assume further that f n is obtained from theorem 6.16. Notice that 
in employing theorem 6.16, v and % are needed. We take v n as v, and let % be determined from r) n 
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in the inductive process, as shown below. Let Vh n and w n be the gradient and divergence free parts 
ofu n , and let z n = w n o (id+V/ n ) o j3 n . Finally, assume that 

II z n ||sj II w n 11 5 5 1| F7h n || s , || u n || s , || /3 n || s , || /3 n || s , || v n || s < R (7.10) 

and 

|| Zn 11s—1 5 || VJ n || s _3, || Vh n || s _|, || V/ n || s , || U n || s _3, || f3 n || g _2) || Vn || s _|— Rj (7.11) 

where (7.10) and (7.11) mean that each one of those quantities is bounded by R. 


Remark 7.1. The reason for introducing u n is that ?y n may not be volume preserving, so that u n 
given by f) n = u n o g n may not be divergence-free. We need, however, a velocity that is divergence 
free to get the correct regularity for the pressure, since it will involve div(V lt u) and we need this to 
depend on at most one derivative of u. 


We shall show that if T and the constants R s , R, etc. are correctly chosen, and k is large, one can 
construct rj n +i and the corresponding quantities satisfying the above conditions. This will give the 
desired sequence. For simplicity, we divide the procedure into several steps. 


Step (nl). We begin by noting that has a smooth normal bundle inside H s (Ll, R 3 ), with smooth 

exponential map, both with respect to the L 2 -metric [26]. Thus, given the curve of embeddings r] n , 
if R s is sufficiently small the exponential map gives 

On = (id +Vg n ) o 7 n , 

with V g n € H s (Ll) and 7 „ € P^(D) having the same regularity properties as r] n , i.e., V g n , 0 n £ 
C 1 ([0, T], H s (Ll)) n C 2 ([0, T], H s ~z (11)). In fact, the normal bundle to D^(D) at r/ is 
\ZrjA ~ 1 div ri (H s (Ll,R n ))-, div^ : H s i-a 17 s ” 1 and V^A ” 1 : H i-a H s are both smooth in 77 , as 
shown in [26]. Thus the normal bundle is smooth in the H s topology even though it is normal only 
in the L 2 sense. Therefore, the smooth exponential map on the normal bundle makes sense and it 
is a diffeomorphism in a neighborhood of the zero-section, and from this it follows that V g n and 
do have the stated regularity. It also follows that Vg n (0) = 0 and 7 n ( 0 ) = id. Define /3 n +i = 7 n - 
This also gives v n+ \ by $ n +i = D 1+1 0 /3 n +i- Since /3 n+ i(0) = id, given R s , if T is small we will have 
|| (3 n +1 ~ id || s < R s - Differentiating rj n in time, evaluating at zero, and applying P (see (4.4)), one 
obtains 


4n+l(0) = « n+ i(0) = Pj)n(0) = PU 0 . (7.12) 

By the bound (7.8b) and the continuity in t of v n+ \, we see that if R is taken sufficiently large one has 
II D 1+1 || s < R an d || $ n +1 || s < R- Differentiating r] n twice in time, evaluating at zero, and applying 
P produces (see (4.6)) 

Di+i(0) + P(V„ n+ 1 ( 0 )V n +i( 0 )) = P{ijn( 0 )) - 2 PD Vn+1 (0) V g n ( 0 ). 

In view of (7.8b) again, we can bound S7g n ( 0), and then, using (7.8d) and the continuity in time of 
the quantities involved, we obtain || u n +i || 3 < R and || f5 n +i |L_s< R, provided that R is chosen 

very large. 
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Step (n2). Define po,n+i by solving 4 : 

\ A Po,n+i = - div(V^Un), on (id+V/ n ) o f3 n (Q), 

\po,n+i = 0) on <9(id +V/ n ) o /3 n (£l). 

Since u n is divergence-free, div(Vg n n n ) is in I7 S_1 , and thus it is expected that po,n+i is H s+1 - 
regular. This has to be verified, however, since one is solving a Dirichlet problem on the domain 
9(id+V/ n ) o p n (Q) = 9(id +V/ n )(D), which is not smooth (compare the ensuing argument with 
that involving Hjj in proposition 6.11). Let rj n = (id+V/ n ), and consider the operator which 
acts on functions defined over D (see notation 2.3). If G is defined over D, a computation gives 

A OnG = ((Drj n )- 1 )p(da'yG)((Drj n )- 1 )'p + d^DrkY^dpG, 

where ((Drj n )~ 1 )p are the entries of (Drjn) -1 . Thus, A jj n has the form 

A 9jn G = a^d a0 G + b a d a G. 

Since V/ n € H S+ ^(Q) we conclude that a“^ is in H s+ 2 (D,M 32 ) and b a is in H s ~z (D, R 3 ). Further¬ 
more, since V/ n is small in H s+ 2(D) and s > | + 2, we find that (a“^) is positive definite. We 
conclude that A^ n is an elliptic operator that takes H S+ 1 (Q) into 7P _1 (fl) and, furthermore, that it 
gives rise to an isomorphism A^ n : Hq + 1 ( f2) -+ iLg _1 (D). In particular, the corresponding Dirichlet 
problem for A^ n is uniquely solvable. But finding p 0 ,n+i hr (7.13) is equivalent to solving 

( Ajj n (p 0)Jl +i o rj n ) = - div(V Sn : u n ) o rj n , in D, 

] Po,n+i o rj n = 0, on dQ. 


From the above, we know that this has a unique solution <7o,n+i € H S+1 (Q,), and thus we obtain the 
desired p 0 ,n +1 = go,n+i 0 ( Vn )~ 1 in H s+1 ((id +V/ n ) o /3 n (D)). 

Next, upon differentiating, we obtain 

A Vn%,n+i = — (div(V^ n u ri ) o rj n )' — ((A^)') <fo,n+i> (7.14) 

and a standard computation shows that 


(A^)-=[V g „,%„. (7.15) 

Combining (7.14) and (7.15), and using the ellipticity of A^ n , we obtain g 0 n+1 G H s ~*(n, R)), by 
using the regularity of the quantities on the right hand side of (7.14). In fact, we know more: from 
the induction hypothesis and the above constructions, we conclude that 

5o,n+i G C°([ 0 ,T},H s+ 1 (n,R))nC\[ 0 ,T\,H s -h(n,R)). 

The ellipticity of A^ n also implies that || go,n+i || s +i can be bounded in terms of || u n || s and 
II V/ n || s _j_ 3 , and that || g 0 ,n+i || S _I is bounded in terms of || u n || s , || u n || s 3 , || V/ n || s+ 3 and 

II V/ n || s . Such bounds hold in particular at t = 0, when they are given in terms of the quan¬ 
tities just mentioned, but now evaluated at time zero. Thus, as before, continuity in t gives a 
bound on || go,n+i || s +i and || g 0 n+i ||s+i in terms of R, provided that R is sufficiently large, i.e., 
II § 0 ,n+1 || s +i< R and || go,n+i || S _1< R- Invoking (2.1) produces bounds for || po,n+i ||s+i and 
II Po,n+i ||s+i in terms of R as well. 


4 Note that the subscript zero in po,n+i is used to indicate that po,n+i is the interior pressure and should not be 
confused with the step of the iteration which is n. 
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Step (n3). With u n +i and go,n+i obtained above, we use theorem 6.16 to solve the /-equation with 
initial conditions /( 0 ) = 0 and /( 0 ) = A“ 1 divuo| (9 ^ (recall that / is determined up to constants), 
obtaining f n +\. In doing so, we need to assure that the initial condition /(0) satisfies (6.33). This is 
the case if || Quq || s is taken sufficiently small, as assumed in theorem 1 . 2 . The bounds obtained for 
v n +\ and (fo,n+i determine the constant K 4 in that theorem. Taking R larger if necessary, we can 
take K 4 = R. Note also that J(id+V/ n+ i) = 1 and id+V/ n+ i is an embedding. 

Remark 7.2. Notice that theorem 6.16 does not say anything about the uniqueness of / n +i, thus 
at this point we let f n +i be any of the (possibly more than one) solutions given by that theorem. 
We will eventually show that uniqueness does hold for the desired equation, i.e., (1.1). 


m 

on 


(id+V/ n+ i)(fi), 

3 (id+V/„+i)(n), 


Step (n4). Obtain h n+ \ solving 

f A/i n+ i = 0, 

lifer = < (V ^ +1 + ^n+i v /n+i + v n+l ) o (id +v/„ +1 )-\iv n+1 ) 

where N n+ i is the unit normal to <9(id +V/ n+ i)(fl). This gives V/i n+ 1 G if® ((id +V/ n +i)(fi)) and 

• 3 

V/i n+ 1 G H s ~2 ((id +V/ n+ i)(f2)). We argue as above to conclude that || V/i n+ i || s and || X7h n+ i || s _ 3 

are bounded by R. Here, again, in producing the estimate we use elliptic theory, so some of the con¬ 
stants involved depend on the domain, namely, on (id+V/ n+ i)(ll), but these constants are controlled 
as before. 


Step (n5). Set rj n+l = (id+V/ n+ i) o f3 n+ i. By construction this is a volume-preserving embedding, 
so the velocity u n + 1 given by rj n+l = u n +i °V n +i is divergence-free and has the same regularity and 
bounds as u n . However rj n+ i is not yet the embedding we are seeking to conclude the n-th step. 
The reason to introduce it is to obtain u n+ 1 , since it will enter in the equation for z below, but it 
has to be defined on (id+V/ n _|_i) o f3 n+ i(Q) (on which u n is not defined). Using as input V/i n+ 1 , 
u n+ 1 , V/ n+ 1 , and /3 n + 1 , and the initial condition Puq, consider the following equation for z n+ \ (or, 
equivalently, for w n+ \\ compare with (7.4) and (7.5)): 


Z n +1 — Qr? J1+1 ((V Un+1 ) r?ri+1 (z n+ 1)) Pv n +i((y Zn+iOTj-l^Vn+l 

+ (V/i n+ 1 O rj n+1 )) + VH n+ 1 o rj n+1 , 


(7.16) 


with initial condition z(0) = Puq, where uq is the initial velocity given in the statement of theorem 
1.2, and VH n+ \ is divergence-free and has normal component equal to (w n+ 1 , (Nn+i°Vn+i)'°Vn+i) + 

(V«„ +1 w n +i) N n +i). 

We shall first show that (7.16) has a solution in i7 s_1 , and then that z n +i in indeed in H s . 

Although the boundary oirj n+1 (Q) is not smooth, it is sufficiently regular to guarantee that the 
operators and Qrj n+1 are bounded on H s (rj n+1 (fl)). In fact, V/ n +i G H s+ 2 (H) and /3 n+ i G 

V s ^{yt). Since /3 n+ i(dQ) = dft, rj n+1 (dQ) = (id+V/ n +i)(9fi). Thus, because V / n+ 1 G H S+ %(Q), 
dfj n+ i(Q) can be written locally as the graph of an H s+1 (dft) function and, therefore, the normal to 
drin+i (U) is P s -regular. 

The term PV U)n+ 1 V/i n +1 is in H if w n . |_i is in H s ~ l because V/i is in H s , and P is an operator 
of order zero. The term Q'Vy m+l vj n+ \ is in H s ~ l if w n+ i is in H s ~ l and Pw n+ 1 = w n+ i. Indeed, if 
w n + 1 is in the image of P we can write 


Q^u n +I w n+1 — \Q 1^7 U n j r \\ w n+ 1- 
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The commutator [Q, Vu n+1 ] is a zeroth order operator depending on first derivatives of u n +i ; since 
u n+ 1 is IP-regular, we obtain that u n+ 1 w n+ \ is in IP -1 . 

Therefore, (7.16) can be viewed as a ODE for z n + i in Pfj n+ 1 {PH s ~ 1 {rj n+ 1 {^l))), i.e., the space 
of IP ^ 1 vector fields over D of the form X = W o rj n+1 with PW = W. Indeed, an element 
X € Prj n+ 1 (PH s ^ 1 (T] n+ 1 (Q))) is of the form X = P- ]n+i Y , with Y £ PiP _ 1 ( 77 n+ 1 (fi)). But then 
y = IT or/ n+1 , with PIT = IT, thus 

* = ^ n+1 D = {P(Y o rj-l,)) o fj n+1 = {PW) o rj n+1 . 

The right hand side of (7.16) depends linearly on w n +1 , and because composition on the right is a 
smooth map (see, e.g., [22, 26, 46]), we conclude that this ODE has a IP -1 solution z n+ \ for a small 
time interval T. Notice that i n +i is also in IP -1 . 

Letting w n +1 = z n +i °rj ~ +1 , it follows that u; n +i is in IP -1 , that Pu; n +i = u; n +i, and that u; n +i 
satisfies 
du] 

Wn+I + p ( v _ t+iWw+ 1 + V Wn+ 1 V/t n+ i) + Vff n+ 1 O Tj n+1 = 0 in 1J {t} x r) n+ 1 (D). (7.17) 

0<t<T 

We now show that is in H s . In order to do so, suppose first that V/i n+ i, u n + i, Vn+h and 

uq belong to H N+l , where N is some big number larger than s. The above ODE argument then 
produces z n+ i and w n+ \ in H N . We shall establish the following a priori bound 

II W n+ 1 || s < C(|| W n+ i ||o +( 1 + || Vn+l III) II u 0 ]|s) 

X eCB+lldill^/o^PlI^n+ilHKIK+ilU+llv^n+ilU). ^ 7 ' 18 ^ 


Once (7.18) is established, one can take a sequence {'Vh n+ i j j,u n+ i t j,rj n+ 1 j,uo ! j } < ^ =1 of H N+1 func¬ 
tions converging in iP to the original V/i n + 1 , u n + i, V n +i ; u o, and then (7.18) implies that the 
corresponding H N solutions w n -|_ij converge in IP -1 to an element that is in fact in IP; this 

will indeed be a solution of (7.17) because we already know that is has a solution in IP -1 . 
Let y = curlu) n+ i. Taking the curl of (7.17) and using the fact that curlV = 0 gives 

|| + Vu„ +1 y+ [curl, Xu n+1 ]w n+1 + [curl, V Wn+ 1 ]Vh n+1 = 0 in |J {t} x fj n+ 1 {n). (7.19) 

0 <t<T 


Here we used that fact that for any vector field X, curl PX = curlX since curlQ = 0. Composing 
(7.19) with rj n+l leads to 

(y°Vn+ 1 )'= —([curl, Vu„ +1 ]wiji+l) °Vn+l — ([curl, V w „ +1 ]V/ln+i) O Tj n+1 . 


Thus 


[curl, V Sn+1 ] 
derive 


y o rj n+ i = {y o 7/ n+ i)(0) - / ([curl, V nri+1 ]u; n+ i) o rj n+1 

Jo 

- / ([curl, V %+1 ]V/i n+ i) ofj n+1 . 


(7.20) 


is a first order operator depending on first derivatives of u n+ \. Thus, using (2.1) we 


|| ([curl, Vu n+1 ]u; n +i) o y n+l || s _i< C || u n+ 1 ||,|| w n+1 || s (1+ || y n+1 ||*), (7.21) 

Similarly, [curl, V Wn+1 ] is a first order operator depending on derivatives of w n+ 1 , hence 
II ([curl, Viu n+1 ] V/l n _|_i) O T ] n +1 1 1 S' — 1 — C || IUII IUn +1 11 s (1+ || Vn+l ID- 


(7.22) 
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Combining (7.20), (7.21), and (7.22) produces 

,, (7-23) 

where we used that (y orj n+ 1 )(0) = y{ 0) = curlito- But 

II CUllw n+1 I|s —1 = 11 y IU—1=11 y°Vn+l°Vnll IU—1< C ( l + II Vnh IIS) II y°Vn +1 IU- 1 , 
which combined with (7.23) gives 

|| curl w n +\ || s _i 5s C*(l+ || f] n _ |_j || s ) || uq || s 

ft (7.24) 

+ ^(1+ II TJ n -1-1 L) / (Id - II Vn+l L)(ll u n+1 II, + II Vh n+1 || s ) || W n +1 || s . 

Jo 

We now use (2.4) to estimate w n+ \ in H s , noting that the terms divu; n+ i and (w n+ \,N), where N 
is the normal to dr} n+ i(£l), do not contribute because Pw n+ \ = w n+ \. We also note that we are 
allowed to invoke (2.4) because rj n+1 is in H s and s >| + 2. Thus, (2.4) and (7.24) give 

II w n+1 ||s — C || W n -|_i 11o TC( 1+ || Us) || Uo ||, 

T || T] n +\ || s ) [ (1+ II Vn+l lls)(ll u n+1 II, + II V/ljj- 1-1 ||s) II W n + 1 || s . 

Jo 

Iterating this inequality now produces (7.18). Note that w n+ \ € H s also gives z n+ \ in H s . As before, 
one gets a bound on || z n+ \ || s and || z n+ \ || s _i in terms of R. This finishes step (n5). 

Now that we have z n+ \ (or w n+ 1 ), V/i n+ i, V/ n +i, and /3 n +i, we define 

r]n+i = id + / {z n+l + Vh n+1 o (id+V/ n+ i) o p n+1 ). (7.25) 

Jo 

Given R s , if T is small, we get || ^ n +i — id || s < R s ■ r/n+i will be in Emb s (0) if it is sufficiently close 
to the identity, and we choose R s accordingly. 

Differentiating r] n+ i we get 

Vn +1 = Zn+l + v/i n+ 1 O (id+V/ n+ i) O /3 n+ i 

= (w n+1 + Vh n+ 1 ) O (id+V/ n+ i) O p n+ i, 

so we let 

U n +1 = t+n+l + Vh n+1 . 

Then u n+ \ is divergence-free and has the same regularity as u n . From the previous bounds, we 
obtain the desired estimates of u n+ \ and u n in terms of R, possibly after increasing R. From the 
initial condition for the /-equation in step (n3), and (7.12), it follows that f] n . (0) = no- 

Differentiating rj n +1 in time twice and evaluating at t = 0, we can control rj n +1 (0) in terms of the 
other quantities of the n + 1st iteration at time zero, which are inductively bounded by a function of 
R s , uo, T and R° s in view of (7.8d). Relabeling the constants and choosing R° s and R appropriately, 
we conclude that the desired inductive bounds hold for the quantities of the n + 1st iteration for the 
desired time interval. 

Finally, a careful analysis of the above steps reveals that, after T and the constants R s , R® and R 
are suitably chosen in order to assure that the quantities of the n + 1st iteration satisfy the inductive 
assumptions, these constants can be chosen uniformly, i.e., independent of n. In particular there 
exists a T > 0 that works for all n, and we get a well-defined sequence {r] n } of embeddings, provided 
that the procedure holds at n = 0, i.e., to start the above iteration, i ]q is needed. Let C be a solution 


V ° Vn+i IU— l < C \\u o ||, +C 


A+ 


Vn -\-1 


^n+1 ||s || ^^n+1 |U) II ^n+l 
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to the Euler equations in the fixed domain 17 with initial conditions £(0) = id, £(0) = Puq, where 
uq is the given initial condition for the free-boundary Euler equations (1.1). For any R s that we 
choose, we can pick T sufficiently small such that || C — id || s < R s . Letting i? be given by ( = i? o 
standard energy estimates for the Euler equations produce bounds for || i? || s and || i? || s _i in terms 
of a constant that depends on T and uq. Hence, we can find a constant Cq(T,uq, R s ) depending on 
T, uq, and R s such that 

II C Id lls^r Rsi 

|| || s < C(T, uq, R a ), || C || s < C(T, u 0 , R s ), (7.26) 

II $ |U-1< C(T, Uq, R s ), || C ||s-l< C(T, Uq, R s ). 

Set rjo = fa = £, uq = i?, V/o = zq = wq = V/io = 0. Minor adjustments have to be made at the 
first step, since /o here is not obtained from theorem 6.16, and rj o(0) = Puq rather than uq, but it is 
clear that this does not hinder the construction of the sequence {rj n }. 


7.3. Convergence of the approximating sequences. Denote by W k, °°([0,T], H s (Q)) the usual 
Sobolev space of H s (£i)-v. alued functions on [0, T] whose derivatives up to order k are essentially 
bounded with respect to the H s ( D) topology. W k,2 ([0, T], H S (Q)) is similarly defined using the L 2 
inner product respect to t € [0, T], 

We start by establishing some further bounds. In the arguments below, the particular form of 
some of the expressions involved will be omitted for the sake of simplicity, since such expressions are 
cumbersome. The relevant information will be the derivative counting. 

Differentiating (6.1) in time, invoking (2.2), recalling that s > | + 2, and using our bounds on /, 
v, and §o, we find that {V f n } is bounded in 77 S_ 3 (D) (with a bound depending on k). 

From step (n4) of the inductive construction, the function h n (t ) satisfies 


^ h n -\-1 0, 


( dhji+i 

8N n+1 


((V/n+1 + D Vn+1 V f n + 1 + V n+ i) o (id+V/ n +i) \N n+1 ) 


in U {t} x (id+V/ n+ i)(D), 

0 <t<T 

on |J {t} x 3(id+V/ n+ i)(D), 

o <t<T 


Differentiating twice, we see that we can bound || V/i n |J s _ 3 in terms of: || V / |j s _ 3 , (which was just 

estimated) together with a constant depending on the domain (id+V/ n _)_i)(D), which was handled 
as in section 7.2; and other quantities that have already been bounded. We obtain therefore an H s ~ 3 
bound for the sequence {Vh n }. 

A similar argument using the equation for z n in step (n5) of the inductive construction shows that 
we can bound {z n } in H s ~ 3 . 

We now establish the convergence. It will be implicit in the arguments that we will be seeking 
convergence of some sub-sequence. 


Convergence of {ij n }: From (7.8a), (7.8c) and (7.8e), we find that the sequence {r] n } is bounded in 
VE 2 ’ 2 ([0, T], H s ~2 (D)), and thus it has a subsequence, still denoted {r] n }, converging weakly to a 
limit r/. Also, r] n and r/ n are bounded in L°°([0, T], LP(D)), thus they have a weakly convergence 
subsequence. We conclude that rj € W 1 ,oo ([0, T], H s (£l)). In particular, rj € C°([0, T], TP(fl)). A 
similar argument yields ij € L°°([0, T], H s ~ 2 (fl)), from which we find that rj € H / 1 ,oo ([0, T], H S (Q)) n 
W 2 ’°°([0,T], H s ~ 1(D)). 

Furthermore, || ij n (t) — rj n (t') || s < R\t — t'\ in view of (7.8c), and {rj n } has compact closure in 
i7 s_ 1 (D) because of the the compactness of the embedding H s C i7 s_1 . Hence, by the Arzela-Ascoli 
theorem, the convergence r) n —>• rj occurs in C 0 ([0, T], i7 s_ 1 (D)). Now, since {rjn} is bounded in H s , 
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interpolating between H s and H s ~ l shows that r] n —>• rj in C'°([0, T], H s ~ s (£})), where 5 > 0 is some 
fixed small number. A similar argument using (7.8e) gives that i) n —> rj in C°([0, T\, After 

interpolation, we have in fact rj n —>• r) in C°([0, T], H S ~ S (Q)). Therefore, rj n —>• r] inC 1 ([0,T],i7 s -' 5 (fi))) 
Next, from the definition of rj n , the inductive bounds and the bounds established above on z n , 
V h n , and V/, we find that { V n } is bounded in H s ~ 3 . Thus, invoking once more the Arzela-Ascoli 
theorem, we find that rj n converges in C°([0, T], H s ~ 3 ~ s (£1)), where 6 > 0 is a fixed small number. 
By interpolation and (7.8e), rj n converges also in in C°([0, T], 

Summarizing: 

• r/ n -> rj in C rl ([0, T], H s ~ 6 (n)) n C7 2 ([0, T], H s ~^~ s (n)), and 

rj G B 7 h°°([ 0,T},H s (n)) n W 2 ’°°([0, T], H s ~i (0)). 

Convergence of {(3 n }: We can repeat the same argument with the sequence {/ 3 n } and conclude that it 
has a limit in VF 1 ,oo ([0, T], H S (Q)) n W 2 ’°°([ 0, T], H s ~^{£!)). Because is a closed submanifold 

of H S (Q), we have in fact P € W 1 ,oo ([0, T],P* (fi)) n VF 2 ’°°([0, T], TF"! (12)). Because of (7.10) and 
(7.11), we can as before use the Arzela-Ascoli theorem and interpolation inequalities to conclude 
that 

• Pn —>• P in C 1 ([ 0 , T], and 

P G 0,r],2?®(fi)) n W 2 ’°°([0,T],i7 s -i(Q)). 

Convergence o/{V/ n }: Recalling that bounds on || / n || Si g translate into bounds on || f n || . i (see 
estimate (6.3)), we apply an analogous argument to the sequence {V/ n }. Using the bounds and regu¬ 
larity given by theorem 6.16, the sequence is bounded in VU 2,2 ([0, T], H s ~z (0)). Also {V/ n }, {V/ n }, 
and {V/ n } are bounded in L°°([0, T],H s+ % (12)), L°°([0, T], if® (12)), and L°° ([0,T],(12)) , re¬ 
spectively. Therefore, V/ n converges weakly in 1U 2 ’ 2 ([0, T], H s ~2 ( 11 )) to a limit V/ which is in 
L°°([0,T],iF + §(12)) n VW 1 ’ oo ([0,r], i^T s (0)) n W 2 ’°°([0,T],i7 s -t(12)). Furthermore, in light of the 
previously obtained bound on V/ n , we have || V/( 2 ) — V/( 2 ') || s _ 3 < C |2 — 2 '|. Hence, as before, 
combining the Arzela-Ascoli theorem with interpolation inequalities gives 

• V/„ ^ V/ in C'°([0,T],i7 s+ |- ,s (O))nC 1 ([0,T],i7 s - < 5 (ll))nC 2 ([0,T],i7 s -|- 5 (Q)) and 

V/ G L°°([0, T],H S+ l(O)) n W 1 - 00 ^, T], FT S (12)) n W 2 ’°°([ 0 , T\,H S ~ 1(17)). 

The other quantities appearing in (7.25) are handled in a similar fashion. We obtain: 

• u n o (id +V/n) o p n —>• u o (id +V/) o p in C°([0, T], 7F” 5 (12)) and 

uo (id +V/) o p G L°°([0, T], 77 s (12)) n kF 1 ’ oo ([0, T],H S ~^ (11)). 

• go,n ->• go in C'°([0,r],77 s+ 1 -' 5 (12)) and 

go G L oo ([0,T],i7 s+ 1 (O))nW 1 ’ oo ([0,T],i7 s -^(O)). 

• V/in O (id +V/n) O /3n -)• V/i o (id +V/) o P in C° ([0, T ], H s ~ 5 (12)) and 

Vh o (id+V/) o p G L oo ([0,T],77 s (Q))nW 1 ’ oo ([0,T],77 s -t(O)). 

• w n o (id +V/n) o P n — y Vw o (id +V/) o p in C°([0, T], H s ~ & {Ct)) and 

wo (id+V/) o p G L°°([0,T],77 S (12)) n lF 1 ,oo ([0, T], 77 s_ 1 (12)). 
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Notice that all the limit quantities satisfy the same bounds as the corresponding sequences. In 
particular, V/ satisfies the bounds given by theorem 6.16. 

7.4. Solution. With the above information, we can pass to the limit in (7.25) obtaining 

77 = id + f (w + V/i) o (id+V/) o (3. 

Jo 

77 is volume-preserving and its velocity u given by f] = u o 77 agrees with u. Also, 77 necessarily 
has the form 77 = (id+V/) o /? and / and f3 have the above regularity properties. In particular 
77 £ <^f(n). Moreover, in light of the way w and Vh were constructed (see section 7.1), and the 
previously established convergences, u satisfies 

da 

P(— + V u u) = 0. (7.27) 

We also know that V/satisfies (4.12a), withp satisfyingpo(id +V/) = poo(id+V/)+Af/o(id+V/), 
where po o (id+V/) = qo (compare with (4.2) and (4.3)), so in particular p = A on dp(Q). From 
the above convergence, we know that po is in H s+1 (Q(t)) and A, being third order in f\ a Q, is in 
H s ~ 1 (dQ(t)), so that Ah £ 7 P _ 2 (f 2 (t)), and hence p £ H s ~z (12(f)). 

It remains to show that (1.1) is satisfied, which is not immediately obvious since we did not solve 
(7.2) in the iteration, but rather P of that equation, i.e., (7.27). We will now show that Q of (7.2) 
follows from the equation for / that we solved. 

Equation (7.27) gives + V M u = + V«u), so there exists a function x such that 

d-u 

— + V M u = V%- 

We need to show that x = ~Pi where p is as in (4.3). As ij o rj^ 1 = jjf + V u u, we have 

rj = 'Vx°V- (7.28) 

On the other hand (7.6) also holds, and because it is a fixed point of the above iteration, / in this 
case does satisfy (4.12a). Therefore, differentiating (7.6) in time twice, decomposing according to 
(4.11) and using (4.12a) and (7.28) gives (compare with (4.6)): 

<2(id -LL^P)(y X o rf) = Q(id -LL^P){-Vp o rj), 

or 

Q(id -D 2 fL^P)(V x o rj) = Q( id -D 2 fL^P)(-Vp o rj ), 
since QP = 0. We can write this as 

{Qrj -1 - (D 2 fL^P)^) (Vx) = (Qr 1 - (D 2 fL^P)^) (—Vp). 

Since Q is the identity on its image, for / small the operator 

Qri -1 - (D 2 fL^P)^ 

is invertible on Q(H s (rj(Q,))), and therefore Vx = —Vp, as desired. 

From the regularity of / and we have Ah 0 (id+V/) o /3 £ L°°([0, T], H s ~ 2 (11)), and hence 
p€L°°([0,T],^-5(Q(t))). As p\ am = A, we conclude that there exists 

77 £ VF 1 ,oo ([0, T], H S (Q)) n VF 2 ,oo ([0, T], H s ~^(Q)) 

satisfying (1.1), as desired. Uniqueness of 77 now follows from the uniqueness of the decomposition 
77 = (id+V/) o f3 given by the exponential map near the identity. This implies uniqueness of p in 
view of (4.3c), (4.3d), (4.3e), and (4.3f). 
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To finish the existence part of theorem 1.2, we point out that 77 € W 1 ,oo ([0, T], H s (£l)) implies 
77 € C 0 ([0, T], H S (Q)), and thus 77 has the stated regularity. 

7.5. Proof of theorem 1.2: existence. Define a new time variable by t = ar, where a > 0 is a 
constant that will be chosen. Define r] a by rj a (r) = 77 (f), i.e., 77 a (r) = rj^ar). Then f/ a (r) = a 2 fj(t) = 
a 2 fy(ar). Using the equation for 77 , i.e., (1.1a), we have f/ a (r) = —a 2 Vp(t) 077 (f) = —a 2 Vp(ar) or](aT). 
So if we define p a (r) = p(t ), i.e., p a (r) = p(ar), we obtain 77 a (r) = —a 2 Vp a (r) o 77 a (r). Then letting 
7r a (r) = a 2 p a (r), we finally obtain 

Va(r) = -V7r a (r) o 770 (r). (7.29) 

Multiplying (1.1c) by a 2 gives a 2 p(t ) = a 2 p(ar) = a 2 p a (r) = 7 r a (r) = a 2 uA{t) = a 2 KA{ar). Thus, if 
we define M a (i") = M(f) = M(ar) and = a 2 u, we have 

vr a (r) = K a M a (r) on 77 a (r)(< 9 D). (7.30) 

Equations (7.29) and (7.30) are of the form (1.1) (with u a (r) defined accordingly) with a coefficient 
of surface tension given by n a . If k > 0 is fixed, not necessarily large, we can then choose a 2 
large enough so that K a is sufficiently large as to apply the result of section 7.4, provided the other 
assumptions can also be accommodated. This is discussed below. We therefore obtain a solution 
(r ] a , 7 r a ). Reverting back to the original variables, this yields a solution to the original problem for 77 
and p with a given k > 0 . 

The result in section 7.4 assumes that dQ has constant mean curvature. We shall show that if we 
are interested only in part ( 1 ) of theorem 1 . 2 , this assumption can be removed as well. 

First, without such an assumption, we do not necessarily have Aqq > 0 so the proof of lemma 6.4 
has to be altererd. Before we used Adn > 0 to show the positivity and invertibility of — A d u — ^AdnA. 
We then used this result in proposition 6.11 to construct an evolution operator associated with A K (t) 
(see equation (6.11)). In the present case, we consider the operator —Ad u instead of —Ad v — \AqqA. 
As discussed in lemma lemma 6.4, —Ad v is positive and invertible, so we still obtain the operator S 
(see (6.10)). This of course gives and extra term in the operator A K (t), namely, ^AqqAS^ 1 . But 
this will be a bounded operator and therefore we still obtain an evolution operator from theorem 2.5 
(see the last statement of theorem 2.4). 

Second, when the mean curvature of D, Man, is not constant, then the equation for / will contain 
the additional term ftMan (see equations (6.1), (5.19) and remark 6.1). This extra term is simply an 
extra inhomogeneous term that can be absorbed into Q. (see equation ( 6 . 11 )). 

The result in section 7.4 also assumes V/(0) to be small, i.e., the gradient part of the initial 
velocity, Quo, has to be small (we do not have to worry about V/(0) being small since V/(0) = 0). 
We now show how this assumption can be removed. 

The assumption that V/(0) is small is used in step (n3) of section 7.2 to guarantee that /(0) = 
A ” 1 divuo | ar2 is small. That /(0) is small is used in theorem 6.16 (see (6.33)) in order to obtain 
estimate (6.36). However (6.36) still holds if /(0) is not small (we explain below why the ensuing 
argument was not used in the proof of theorem 6.16). I.e., in theorem 6.16, assume that instead of 
(6.33) we have 

II fi \\ 8 + i td < K 3 . (7.31) 

In what follows, we continue to assume that k is large since, as showed above, the problem for 
arbitrary k > 0 can be reduced to that of large k via a rescaling. 

As in the proof of theorem 6.16, we invoke (6.35). The estimate of the term lA{t,T)Q(r ) does 

not rely on (6.33), so this term yields sup 0<r<T || G(t) || s+ i q as before. 
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For the term U(t, 0)z(0), first notice that (6.36) corresponds to the first component of z, i.e., z\. 
Write 


W(t,0)z(0) 


'Un(t,0) Ui 2 (t,0)\ fo\ 

w 2 i(t,0) u 22 (t,o)J V/J’ 


where z(0) = (0,/i). The first component of the above is Ui 2 (t, 0)fi- Recall that U( 0,0) = I, so 

S+ 1 

£R2(0,0) = 0, and also U(t,r) is strongly continuous into Y = H 0 2 (dQ) (see theorem 2.5 and 
section 6). Thus, with f\ and k given, we can choose T (and hence t) so small that 


Wl2M)/i 



C(Kq)K 3 

K 


Therefore, estimate (6.36) still holds without the assumption that V/(0) is small. 

The above argument was not used in the proof of theorem 6.16 because it produces a time interval 
[0, T k ) that shrinks to zero as k —> 00 , so the corresponding existence result and estimates would not 
apply to the limit n —>• 00 . 

The other part in the proof of theorem 6.16 where (6.33) has been employed was (6.37). It is 
clear, however, that the argument following (6.37) still holds if (7.31) replaces (6.33). Indeed, under 
(7.31), (6.37) becomes 


h || s +i a < C(K 0 ) || h 


<C(K 0 )K 3 + C(K 0 )T sup 

0<t<T 


5 +i a +C(K 0 )T sup 

2 ’ 0 <t<T 

&(t) || 


Q{t) 


s +f>0’ 




(7.32) 


and we can still choose £ large enough so that the right hand side of (7.32) is less than £. 

Finally, (6.37) was also invoked when we derived estimate (6.46). But again, the ^ factor is not 
needed here since we only need the right hand side of (6.37) to be bounded in order to obtain (6.46), 
which is the case in light of (7.32) (see the paragraph after (6.46)). 

Inspection in the proof leading to the existence part in section 7.4 shows that the remaining 
arguments are the same without the assumption that Quq is small. This establishes the proof of 
theorem 1.2. 


Remark 7.3. Notice that these arguments are consistent, in the following sense. We obtain a 
solution that exists for 0 < r < T, or, in the t variable, 0 < -t < T, i.e, 0 < t < aT■ So, if we 
take the limit a —>• 0, so that K a —» 0, the interval of existence shrinks to zero, as it should since 
the problem is not well posed when n a = 0 and the Taylor sign condition, which we do not assume, 
does not hold (it turns out that T also depends on a, so this idea of consistency with a —>• 0 is more 
complicated than just stated, but on a heuristic level we see that we obtain what is expected). 


8. Proof of theorem 1.2: convergence 

Here we establish the convergence part of theorem 1.2, and thus we assume the corresponding 
hypotheses and notations throughout. Some of arguments below resemble those of theorem 5.1 in 
[25] and theorem 5.5 in [23]. From now on it is convenient to re-instate the subscript k. 

Let [0, T k ) be the maximal interval of existence for the solution (r] K ,p K ) found above. Let T K = 
min{T, T K }, where we recall that [0, T] is an interval on which the solution to (1.2) is defined. We 
henceforth consider the quantities r] K and £ on [0,77c)- Assume also that T is chosen such that (7.26) 
holds. 

As in the calculations leading to (4.6), we differentiate r) K twice in time to obtain 
f)K-h = D 2 fn o + (V/ K + 2 D VK Vf K + D 2 VKV Vf K ) O P K , 


( 8 . 1 ) 
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where we used that ft K = v K o (3 K , so that 


, dv K 


Pk = (-QjT+Vv K V K ) o/ 3 k 


Integrating (8.1): 


i)k - Pk = - Pu 0k + 


/W 

Jo 


+ 2 D VK Vf K + Dl v Vf, 


k) ° Pk + f 
Jo 


D 2 f K o /3 k / 3 k , 


where we used that rj K ( 0) = uq and P K { 0) = Puq. Write (8.3) as 

Vk Pk — Qu Ok d - Rki 

where 


Rk. = 


fr„. 

Jo 


with 


Tk = (V/ K + 2 D v yf K + D 2 X7f K ) O f3 K + D 2 f K o fjJ K 


( 8 . 2 ) 

(8.3) 

(8.4) 

(8.5) 

( 8 . 6 ) 


On the other hand, (4.4) gives 

V k-Pk = (V/ k + D VK Vf K ) O p K , 
so that the estimates of section 7.2 give 

II v k - Pk lls <|| v/ K o /3 k || s + || {D v y f K ) o p K || s 

< c || vA II. (i+ II Pk IID + c || v K || s || v/ K |U+i (i+ || & ||2) 

< eg) + C^) < C(R) 


(8.7) 


-y/ft K s/K 

where C(R) is a constant depending on R, and R is as in section 7.2. Combining (8.7) with (8.4) 
and our assumptions on Quo K leads to 

C(R) 


Rk ,< 


\[P 


( 8 . 8 ) 


Remark 8.1. It is important to notice that (8.8) follows from the specific relation between R K and 
Vk — Pk j be., equation (8.4). In fact, we could not try to estimate R K term by term from (8.5) and 
(8.6) since some of such terms, e.g. V/ K , may not even be in H s . 

Equations (8.4), (8.5) and (8.6) also give 

Vk ~ Pk = r K . (8.9) 

Recall that we denote by ( the solution to (1.2). From (8.4) we have 

C - Vk = C - Pk ~ r K , 

so that 


C - Pk = $0 - UOk + [ (C - Pk) ~ R, 
Jo 


( 8 . 10 ) 


after recalling (8.5). 

It is well-known (see, e.g., [26]) that (1.2) can be written 

C = Q(V#0) o c, 
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and, since i9 = £ ° C U it follows that 

c = (Q(v CoC _ 1 cor 1 ))oc ( 8 . 11 ) 

On the other hand, from (8.2) we find, using that Qv K = 0, that 

Q(^o/ 3 ~ 1 ) = Q(V Vk v k ), 

and thus 

+ V Vk v k = /3 k o p- 1 = Q(/3 k O p- 1 ) + P(P K O p- 1 ) = Q(V Vk v k ) + P(p K o Pp 1 ). 
Composing with /3 K and recalling (8.2) once more: 

k = (Q(V„ k u k )) o P K + (P(P K o p- 1 )) o p K , 


so (using notation 2.3) 

Ac = (Q(Vp K oK'P« ° fc 1 )) ° & + PpJ- 

Using (8.11) and (8.12) into (8.10) produces 


C - i] K = i?o - u 0k - R k - / Pf j J3 i 


f 


+ 


f ((Q(v^ oC _iCo r 1 )) ° c - (q(v^ d/3 - i/ 3 « o p- 1 )) o p K 

J 0 


(8.12) 


(8.13) 


The term Q(V fi KO/ g-i Ac ° fin 1 )) ° Pk, is in H S (Q) because P K € T>^(fl) and so is (Q(V^_iC° C *)) ° C- 

Integration of (8.12) shows that /[[ Pp K 'k is in H S (Q) in that Ac € H S (Q). We can, therefore, estimate 
(8.13) in H s and we proceed to do so. 

Denote by the tangent bundle of T>*(fl) and let Z be the map 

Z : TT>® (fi) -> H s (n, R 3 ), 

Z(Z,X) = (Q(V Xo t- ilof^oe 


In [26] it is shown that Z is a smooth map. Since the image of (CC) is compact, Z is uniformly 
Lipschitz in a neighborhood of (CO- Thus 

ii z((, o - z(p K j K ) || S < cdi c - a, iu + ii c - k y • (8.i4) 

Combining ( 8 . 8 ), (8.13), (8.14), and our assumptions on uq k , it follows that 

II C — 7 k ||s < — 7=^(1 + t)+ || f Pp K Pn IU +C (1 + t) [ || C ~ 7k |U> (8.15) 

V K Jo Jo 

where we used the fact that the term || C, — P K || s in (8.14) can be estimated in terms of || £(0) — 
P K ( 0) || s + || C — P K || s because of the fundamental theorem of calculus. We also used the fact that 
II C - Ac |U<|| C ~ 7 k IU + 7 =- Thus, iterating (8.15), 

II C-7k II, < (^(1 + t)+ || J*PpJ* IU )e C(t+t2) . (8.16) 

Integrating (8.12) and using the estimates of section 7 yields 




52 


DISCONZI AND EBIN 


By the definition of T K , || C IU remains uniformly bounded on [0,7^), and therefore the same holds 
for || f] K || s in light of (8.16) and (8.17). We conclude that [0,7^) is not the maximal interval of 
existence of the solution of (1.1), and therefore T K > T. Moreover, this conclusion holds for all n 
sufficiently large since (8.16) and (8.17) hold for all k sufficiently large. 

.. 3 

Next we proceed to show convergence. First we estimate the term Pp K /3 K in H s ~ 2 . Using (8.2) 
and (4.12c) we get 


m ° fc 1 ) + L^P(2D v „Vf K + D 2 VKV Vf K ) 

+ L^P(L k Q(V Vk v k )) = -L^P(Vp K o (id+V/«)). 


(8.18) 


All terms in (8.18) are in H s 2 (D), and we have 


L^P(2D VK Vf K + Di KV Vf K ) || s _3< 


C(R) 


K 


and 


L^P(L k Q(V Vk v k )) || s _3< 


C{R) 


2 K 

after using PQ = 0. 

Writing q K = p K o (id+V/ K ), we have 

Vg K = Vp K o (id+V/ K )(id+D 2 / K ), 


so 

Thus, 


Vp K o (id +V/«) = Vq K + 0(D 2 f K )Vq K . 
P{Vp K o (id +V/ K )) = P(0(D 2 f K )\7q K ), 


so that 

II P(Vp K o(id+V/ K )) || s _3< 

2 K 

These estimates and (8.18) imply 

II Ppj% || s _|< (8.19) 

Thus (8.16) and (8.19) imply that || ( — || s _3 —> 0 as k —>• 00 . The interpolation inequality can 

now be invoked to conclude that || ( — || s _ 5 —>• 0 for any fixed 5 > 0. We thus also conclude that 

that || JqP(3 k P k IIs—<5—^ 0. 

Let Jrj be a standard Friedrichs modifier in S7. Recall that if h is defined in U, J a h is given by 
(see, e.g., [64]) 


J a h = RJ a Eh, 


( 8 . 20 ) 


where E : H r (Q) —> H r (M. n ) and R : H r (M . n ) —> H r (Q) are, respectively, the extension and restriction 
operators, and J a is the Friedrichs modifier in M n . Notice that, here in contrast to (2.3), one does not 
lose a half-derivative upon restriction since U is of the same dimension as M n . Denoting by || • || s rt> 
the Sobolev norm in M n , standard properties of J a give 


Jrrh 


C „ , 

1 < —3- || n 

(7 2 


( 8 . 21 ) 
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if h is in H s 2 (1C). From (8.20) and (8.21), 


l 


Ja / P/3 k Pk |U =|| RJ*E 


f'PpJn 

Jo 


,< C II J a E 


( Pp K PK ||s,I 

Jo 


C f t 

— T II E PpuPn || s R ™ 

(72 Jo 2 

c r f ■■ ct 

Ppjn || a _|< — || Pp„p K II§—| 

2 a 2 2 


< 3 


/' 

Jo 


Thus, in light of (8.19) 


(72 Jo 


rt CT C(R) 

J(T / P^Pk ||s< 3 “ II Pp K Pn || s _|< 3 —, 
Jo 


(72 


(7 2 X /K 


( 8 . 22 ) 


where T was absorbed into C{R). Consider 


J^ P pJk = J*PpJ«~ J ° I PjJn + Ja I Ppjn- 


f Ppj. + Ja f 

Jo Jo 


and choose a sequence of cr’s depending on k by a = a K = Then a K —>• 0 as k —> oo and (8.22) 
gives 

rt 


l 


J&K / PPkPk II 


C(R) 

1 1 
K4 


But 


and thus we conclude that 


lini II [ Pp K P k J<r [ P/3 k Pk lls — 
Jo Jo 


cr—>0 


linr || f Pp K /3 K || s = 0. 

K^-OO J Q 

Combining this with (8.16) gives || f] K || s —>• 0 as k —> oo. This immediately produces || (—r] K || s —>• 0 
as well since £(i) = £(0) + J^ £, and similarly for r] K . From the regularity of rj K the convergence is as 
stated in theorem 1.2, finishing the proof. 


8.1. Proof of corollary 1.8. From our estimates, we immediately have 

|| V f K 11 s _(_ 3 < , 

^2 K 

and 

II V f K || s < — -j=i 
V't 

so that V/ K —>• 0 in H s+ z and V/ K —>• 0 in H s , as stated. Combined with decomposition (1.7), this 
gives 

II Vk Pk l|s ^ o, 

so that /3 k —>• C in E[ S (Q) in view of the convergence part of theorem 1.2. Similarly, one gets the 
convergence /3 K —>• ( in H S (Q) from (4.4). 
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